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ABSTRACT 


The  Sacramento  Valley  water  quality  model  simulates  mathematically  tne 
monthly  flow-weighted  mean  electrical  conductivity  (EC)  by  reaches  In 
the  Sacramento  River  system.   Using  readily  available  flow  and  EC  data, 
the  rrodel  vjas  developed  arc'-'"-.i  the  system's  three  ^ereril  *-ype2  of  f  1  t.-,' : 
unregulated  ti'lbutary  inflow;  reservoir  impoundments  and  releases;  and 
valley  floor  accretions  or  depletions.  The  latter  was  further  broken 
down  into  storm-water  runoff,  irrigation  return  flow,  municipal  and 
industrial  waste  discharges,  diversions  of  streamflow  for  beneficial  use, 
weir  spills,  and  unmeasured  Inflows  and  losses. 

Plow  and  EC  data  for  all  inflows,  diversions  and  other  losses  were  com- 
bined in  a  series  of  equations  in  which  mean  monthly  EC  was  weighted 
with  quantities  of  flow.   Model  outputs  were  flow-weighted  mean  EC  in 
micromhos,  salt  loading  in  tons  and  corresponding  rates  of  streamflow  in 
acre-feet  at  12  locations  in  the  system  for  each  month  of  a  selected 
time  period. 

The  model  was  verified  by  comparing  computed  EC  values  with  prototype 
values  at  selected  river  stations.   Although  rather  large  deviations 
occurred  between  computed  and  prototype  EC  values  for  some  months,  the 
model  simulated  the  general  pattern  of  EC  fluctuations  at  the  river 
stations.   Excessive  deviations  generally  could  be  attributed  to  specific 
inadequacies  of  prototype  data. 

The  model  was  used  to  estimate  the  probable  future  EC  at  each  station  in 
the  system  under  a  projected  1990  level  of  development.   Results  showed 
an  Increase  in  flow-weighted  average  EC,  but  no  monthly  values  were  great 
enough  to  threaten  any  beneficial  uses  of  the  Sacramento  River. 
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FOREWORD 


Selection  of  the  most  favorable  management  plans 
for  California's  water  resources  relies,  in  part,  on  careful 
evaluation  of  possible  future  changes  in  water  quality.   Such 
changes  may  be  brought  about  by  altered  land  and  water  use 
patterns,  construction  of  dams  and  reservoirs,  or  alternative 
methods  of  reservoir  operation.   Foreseeing  the  probability 
of  adverse  water  quality  conditions  will  give  planners  oppor- 
tunities to  minimize  or  even  eliminate  some  potential  problems. 

A  computer-programmed  mathematical  model,  which  can 
simulate  a  selected  indicator  of  water  quality,  is  a  valuable 
tool  for  estimating  water  quality  under  future  conditions. 
A  model  developed  by  the  Department  of  Water  Resources  as 
part  of  a  special  water  quality  study  simulates  mean  monthly 
electrical  conductivity  in  the  Sacramento  River.   Electrical 
conductivity  was  selected  as  the  water  quality  indicator  for 
this  study  primarily  because  it  measures  total  dissolved  solids, 
or  salinity,  and  because  it  can  be  used  to  estimate  other  min- 
eral constituents. 

The  study  was  conducted  in  I969  and  197O  with 
counties  and  local  valley  irrigation,  reclamation,  and  drain- 
age districts,  all  of  whom  assisted  with  costs  of  computer 
programming  and  operation  and  provided  a  technical  advisory 
committee  to  work  cooperatively  with  the  Department. 


William  R.  Gianelli,  Director 
Department  of  Water  Resources 
The  Resources  Agency 
State  of  California 
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INTRODUCTION 


Systems  analysis  techniques  employing  mathematical  models  program- 
med for  the  digital  computer  provide  a  highly  flexible  means  of 
incorporating  water  quality  considerations  in  water  resource  plan- 
ning.  A  model  permits  prompt  evaluation  of  the  probable  response 
of  a  particular  water  quality  factor  in  the  prototype  system  to  a 
change  in  that  system. 

The  modeling  technique  was  applied  in  a  special  water  quality  study 
of  the  Sacramento  River,  completed  in  July  1970,  that  was  directed 
toward  analysis  of  factors  that  influence  the  mineral  quality  of 
water  in  the  river  system.   Its  overall  purposes  were  to  provide 
the  Department  with  environmental  planning  information;  to  give 
irrigation,  reclamation,  drainage  districts  and  associations,  and 
county  governments  in  the  Sacramento  Valley  the  information  they 
needed  to  evaluate  the  Central  Valley  Regional  Water  Quality  Con- 
trol Board's  proposed  water  quality  control  policy  for  the  Sacra- 
mento River;  and  to  provide  the  Board  the  means  to  develop  an 
adequate  policy. 

The  mathematical  model  developed  in  this  study  was  used  to  analyze 
present  and  probable  future  electrical  conductivity  at  12  stations 
in  the  river  system  between  Keswick  Reservoir  and  Freeport.   The 
model  permits  evaluation  of  the  effects  on  the  river  mineral  qual- 
ity of  changes  in  land  and  water  use.   The  development  and  struc- 
ture of  that  model  and  illustration  of  its  use  for  estimating 
future  river  salinity  are  described  in  this  report. 

Electrical  conductivity  (EC)  was  selected  as  the  water  quality 
parameter  for  modeling  for  several  reasons:   1)  more  historical 
data  were  available  for  EC  than  for  any  other  constituent;  2)  EC 
as  a  conservative  constituent  of  water  is  not  changed  signifi- 
cantly by  time  or  biological  activity;  and  3)  EC  can  be  used  to 
estimate  other  conservative  constituents,  particularly  total  dis- 
solved mineral  concentrations. 

This  study  was  conducted  as  a  part  of  the  Department's  continuing 
role  in  planning  for  protection  of  water  quality  in  water  resource 
management  planning.   The  model  described  in  this  report  is  not  a 
final  product.   It  is  being  continually  updated  and  refined  as  new 
information  becomes  available.   In  future  studies,  other  important 
water  quality  indicators  than  dissolved  minerals  should  be  evalu- 
ated by  mathematical  simulation  of  the  Sacramento  River  system. 
Nonconservative  nutrient  and  organic  constituents,  which  cannot 
be  simulated  by  this  model,  will  be  of  particular  Importance. 
Consideration  of  these  additional  parameters  will  make  possible 
fuller  evaluation  of  the  impact  of  future  economic  development  on 
the  quality  of  water  in  the  Sacramento  River. 


The  Department  will  make  the  model  available  for  farther  use  by- 
other  organizations  that  seek  solutions  to  water  quality  problems 
in  the  Sacramento  Valley.   Information  on  the  machine  program  can 
be  obtained  through  the  Statewide  Planning  Branch  of  the  Division 
of  Resources  Development. 

Area  of  Study 

The  area  studied  (Figure  1)  included  the  entire  Sacramento  Valley 
drainage  area  above  Freeport,  about  24,000  square  miles.   Primary 
attention  was  focused  on  the  valley  floor  where  occur  virtually 
all  agricultural  drainage,  municipal  and  industrial  waste  dis- 
charges, and  other  significant  factors  that  tend  to  degrade  the 
quality  of  water  in  the  Sacramento  River.   This  report  discusses: 

1)  the  Sacramento  Valley  water  quality  model  as  developed  in  19^5; 

2)  the  current  model  as  used  in  this  study;  3)  computational 
methods  and  assumptions  employed;  4)  the  hydrologic  base  period 
considered;  5)  development  of  model  inputs;  6)  verification  of  the 
model;  7)  estimates  of  1990  mineral  quality  of  water  in  the  Sacra- 
mento River  expressed  as  electrical  conductivity;  and  8)  results 
of  the  computations. 

Data  Considerations 

Data  used  to  develop  and  operate  the  model  were  primarily  those 
which  were  readily  available  in  the  Department's  files  and  reports 
or  in  federal  publications.   Agricultural  drainage  data  collected 
under  a  second  part  of  this  study  were  available  for  use  as  guides 
in  assigning  quantity  and  quality  values  to  agricultural  return 
flows  to  the  Sacramento  River. 

As  expected,  available  data  were  not  fully  adequate  for  detailed 
analyses  of  all  important  factors  which  influence  the  mineral 
quality  of  water  in  the  Sacramento  River  system.   Some  hydraulic 
and  other  characteristics  of  the  system  are  virtually  unknown, 
and  their  determinations,  if  they  were  possible,  would  require 
detailed  field  investigations.   For  example,  quantities  of  seep- 
age in  or  out  of  any  given  reach  of  river  have  not  been  determined 
for  any  historic  period.   Yet  such  seepage  can  significantly  affect 
mineral  quality  of  water  in  the  river,  because  seepage  outflow 
removes  salts  and  seepage  inflow  adds  salts. 

Data  also  are  sketchy  on  the  fate  of  minerals  introduced  to  valley 
floor  areas  in  applied  surface  water,  rim  inflows  and  precipita- 
tion. At  certain  times  such  minerals  are  retained  in  the  soils 
and  groundwater  of  the  area  to  be  leached  out  later.   Investiga- 
tion of  this  phenomenon  could  not  be  included  within  the  scope 
of  this  study. 

Notwithstanding  data  inadequacies,  a  usable  mathematical  simula- 
tion of  the  Sacramento  Valley  system  was  possible.   Inputs  were 
provided  for  all  mineral  sources  considered  significant.   Where 
data  were  weak  or  lacking,  value  judgments  were  made. 


SACRAMENTO  VALLEY  WATER  QUALITY  MODEL 


In  1965  the  Department  developed  a  mathematical  model  for  use  in 
estimating  future  mineral  quality  of  water  in  the  Sacramento 
River.   The  model  was  concerned  mainly  with  the  quality  of  water 
in  the  State  Water  Project,  and  hence,  major  emphasis  was  placed 
on  the  quality  of  Sacramento  Valley  outflow.   The  study  reported 
herein,  on  the  other  hand,  is  concerned  with  the  quality  of  water 
at  additional  points  in  the  river  as  related  to  the  previously- 
mentioned  water  quality  control  policy.   It  was  necessary,  there- 
fore, to  completely  restructure  the  model  to  comply  with  the  new 
objectives. 

The  basic  structure  of  the  I965  model-^is  still  valid,  however, 
and  a  brief  description  of  that  model  would  be  helpful  in  under- 
standing the  revised  model  used  in  this  study.   In  both  models, 
mineral  concentration  is  expressed  in  terms  of  electrical  con- 
ductivity (EC).   For  Sacramento  River  streamflow,  total  dissolved 
solids  (in  parts  per  million)  is  approximately  equal  to  EC  (in 
micromhos)  x  70  percent. 

The  1963  Model 

The  1965  model  was  developed  around  three  general  categories  of 
flow:   unregulated  tributary  runoff;  reservoir  impoundments  and 
releases;  and  valley  floor  accretions  and  depletions.   Flow  and 
quality  data  on  all  inflows,  including  natural  flow,  agricultural 
drainage,  and  municipal  and  industrial  waste  discharges,  and  on 
diversions  are  combined  in  a  series  of  equations  in  which  mean 
monthly  EC  is  weighted  with  quantities  of  flow.   Flow  is  expressed 
in  cubic  feet  per  second. 

Unregulated  Tributary  Runoff 

In  deriving  a  method  for  estimating  future  quality  of  unregulated 
tributary  runoff  into  the  river  system,  a  quality-versus-flow 
relationship  was  employed.  A  reasonably  stable  relationship 
exists  between  rate  of  flow  and  mineral  concentration  for  natural, 
unregulated  runoff  from  the  mountainous  watersheds  bordering  the 
Sacramento  Valley.   This  relationship  can  be  expressed  by  a  math- 
ematical equation  in  the  form  C=aQ'^,  where  "C"  represents  elec- 
trical conductivity;  "Q",  the  rate  of  streamflow;  and  "a"  and  "b", 
constants  characteristic  of  individual  watersheds.   Empirical 


1/  Described  in  detail  in  entry  No.  27  of  the  bibliography, 
JTppendix  A. 


C/Q  relationships  were  developed  by  plotting  observed  values  of  EC 
against  corresponding  rates  of  flow  for  the  headwaters  of  the  Sac- 
ramento River  and  its  major  tributaries  at  gaging  stations  located 
where  the  natural  regimen  of  flow  has  not  been  substantially  altered 
by  major  impoundments  or  diversions. 

Reservoir  Impoundments  and  Releases 

The  natural  C/Q  relationship  is  generally  upset  by  impoundment  in 
a  reservoir.   Reservoir  stratification,  incomplete  mixing,  and 
varying  detention  times  all  complicate  the  problem  of  predicting 
mineral  content  of  reservoir  outflow.  An  emperical  approach  based 
on  available  inflow-outflow  records  was  used  to  develop  a  computer- 
oriented  method  for  deriving  a  quality  of  reservoir  releases  from 
inflow  data.   No  general  solution  was  derived  for  this  problem, 
and  each  reservoir  required  individual  treatment.   The  following 
two  basic  approaches  were  employed. 

1)  The  quality  of  reservoir  releases  was  assumed  equal 
to  the  quality  of  water  in  storage,  which  was  deter- 
mined by  successively  mixing  the  flow-weighted 
average  quality  of  inflow  for  each  month  with  the 
average  quality  of  water  in  storage  for  that  month, 
and  then  allowing  for  a  one -month  detention  time  in 
the  reservoir.   The  inflow  quality  was  computed  by 
the  C/Q  equation  and  the  initial  quality  of  water  in 
storage  was  assumed  equal  to  the  quality  of  inflow 
for  the  month  preceding  the  first  month  under 
consideration. 

2)  The  quality  of  reservoir  releases  was  assumed  equal 
to  the  quality  of  water  in  storage,  which  was  deter- 
mined by  the  monthly  running  average  of  inflow 
quality,  weighted  according  to  flow,  over  a  period 
of  time  in  which  the  total  volume  of  outflow  during 
the  period  equaled  or  exceeded  the  volume  of  storage 
at  the  end  of  the  period. 

Neither  of  these  approaches,  by  itself,  gave  acceptable  answers 
for  all  months  and  at  all  reservoirs.   This  was  shown  by  com- 
parison with  recorded  data.   In  most  cases,  better  answers  re- 
sulted from  a  combination  of  the  two  approaches,  often  together 
with  a  factor  to  allow  for  unknown  impoundment  influences. 
Mathematical  explanations  of  the  two  approaches  are  included  in 
this  report  as  Appendix  B,  "Reservoir  Computational  Methods". 

Valley  Floor  Accretions  and  Depletions 

In  addition  to  its  major  measured  tributaries,  the  Sacramento  River 
receives  local  drainage  and  unmeasured  tributary  inflow  as  it 
traverses  the  valley  floor.   At  the  same  time,  various  amounts  of 
water  are  diverted  from  the  river  for  beneficial  use.   The  changes 
in  water  quality  resulting  from  these  factors  are  represented  by 


the  third  part  of  the  model.   Accretions  which  originate  In  the 
valley  floor  account  for  a  large  part  of  these  water  quality- 
changes.   In  the  Initial  study,  historic  water  quality  data  were 
assumed  to  represent  future  mineral  concentration  of  the  agricul- 
tural drainage  portion  of  these  accretions  and  flow  quantities 
to  be  proportional  to  irrigated  acreage  drained.   The  remaining 
unmeasured  accretions  and  depletions  were  accounted  for  by  eval- 
uating the  differences  in  flow  and  water  quality  at  existing 
stream  gaging  stations  along  the  main  stem  of  the  river.   In 
later  studies,  estimates  of  future  accretion  and  depletion  flow 
quantities  were  taken  from  actual  water  supply  studies. 

In  short,  the  I965  water  quality  model  for  the  Sacramento  River 
used  the  foregoing  relationships  to  determine  the  weighted  mean 
electrical  conductivity  of  tributary  inflow  and  valley  floor 
drainage,  as  modified  by  impoundments  and  diversions. 

The  model  was  verified  by  checking  its  ability  to  reproduce  his- 
toric records  at  key  points  in  the  watershed  and  at  the  head  of 
the  Delta.  An  average  difference  of  about  10  percent  was  noted 
between  computed  and  historic  values  of  conductivity.   Wider 
discrepancies  were  found  for  individual  months,  particularly 
when  large  fluctuations  of  flow  occurred  during  the  month.   This 
difference  may  be  related,  at  least  in  part,  to  the  fact  that 
computed  values  were  based  on  average  monthly  flows,  while  the 
actual  Instantaneous  observations  of  water  quality  that  constitute 
the  historic  record  may  or  may  not  represent  average  conditions. 
Nevertheless,  satisfactory  correlations  were  obtained,  though  it 
was  recognized  that  the  model  could  be  improved  by  collection  of 
additional  data  and  refinement  of  empirical  relationships. 

The  model  has  been  used  with  system  operation  studies  to  evaluate 
water  quality  conditions  in  the  Sacramento  River  watershed  under 
projected  1990  development.   Basic  hydrology  representative  of 
both  normal  and  critical  dry  periods  was  used  and  agricultural 
drainage,  waste  discharges,  and  other  accretions  were  estimated 
for  1990  level  of  cultural  development.   Results  demonstrated 
that,  despite  increased  development  and  water  usage,  the  mineral 
concentration  of  Sacramento  River  inflow  to  the  Delta  will  remain 
well  within  established  State  Water  Project  contract  objectives, 
even  during  the  most  severe  dry  period  on  record. 

The  Current  Model 

The  current  model,  used  in  this  study,  may  be  viewed  as  an 
extensive  refinement  of  the  third  part  of  the  I965  model,  the 
part  which  deals  with  valley  floor  accretions  and  depletions. 
Mean  monthly  electrical  conductivity  (EC)  in  mlcromhos  remains 
the  water  quality  parameter.   Flow  in  the  current  model  is 
expressed  in  thousands  of  acre-feet  per  month. 

Since  only  mean  monthly  conditions  were  considered,  certain 
limitations  of  the  model  should  be  noted.   The  model  does  not 


show  short-term  fluctuations  or  instantaneous  values  of  flow  or 
EC  which  would  occur  at  any  time  during  a  particular  month.   Thus, 
it  does  not  indicate  the  dynamic  effects  of  a  large  momentary  slug 
of  municipal  or  industrial  wastes  discharged  into  the  river.   Nor 
does  the  model  account  for  time  of  travel  of  flow  in  the  river  sys- 
tem. But,  because  travel  time  from  Keswick  to  Preeport  is  only 
about  six  days,  this  latter  factor  is  not  important.   The  model 
also  fails  to  account  for  tidal  action  in  the  lower  reach  of  the 
river,  which  would  affect  the  Sacramento  and  Preeport  stations. 

For  the  I965  model,  the  Sacramento  River  was  divided  into  only  two 
reaches  between  Keswick  Dam  and  Preeport.   The  intermediate  station 
was  Hamilton  City.   Nine  Sacramento  River  reaches  below  Keswick 
were  used  for  the  current  model.   Each  of  the  two  major  tributaries, 
the  Feather  and  the  American  Rivers,  was  considered  as  one  reach. 
The  current  model  thereby  could  more  accurately  describe  the  EC  of 
water  in  the  full  Sacramento  River  than  could  the  I965  model. 

The  nine  Sacramento  River  reaches  below  Keswick  Dam  were  designated 
reaches  "A"  through  "l".   The  Feather  River  was  designated  reach  "j"; 
the  American  River,  reach  "k";  and  all  inflow  above  Keswick  Dam, 
reach  "L".   The  system  is  shown  schematically  on  Figure  2.   Reach 
terminals  were  selected  on  the  basis  of:   1)  water  quality  check- 
points mentioned  in  a  proposed  Sacramento  River  water  quality 
control  policy,  2)  additional  checkpoints  required  to  structure 
the  model,  and  3)  availability  of  water  quality  and  streamflow  data. 

Figure  2,  page  8,  shows  the  identification  number  of  each  source  of 
input  to  the  model  and  its  location  by  river  reach.   The  facing  page 
lists  each  inflow  and  outflow  designated. 

Inflow  sources  included  reservoir  releases,  tributary  streams, 
agricultural  drains,  municipal  and  industrial  (M&I)  waste  dis- 
charges, and  any  additional  inflow  and  other  accretions  that 
could  occur  in  each  reach.   The  distinction  between  inputs  des- 
ignated as  "additional  inflow"  and  "other  accretions"  is  explained 
later. 

Outflow  sources  included  diversions  from  the  river,  weir  spills, 
and  any  additional  losses  that  could  occur  in  each  reach. 

Monthly  flow  and  EC  data  on  each  source  were  obtained  either 
directly  from  available  records  or  estimated  where  recorded  data 
were  not  available  or  where  they  did  not  apply.   Flow  and  EC  in- 
puts for  "other  accretions"  could  be  generated  by  the  model. 

A  vital  component  of  the  current  model  was  an  "accretion  sub- 
routine", in  which  the  quantities  and  qualities  of  valley  floor 
accretions  were  computed  for  each  river  reach.   Inputs  were 
monthly  flow  quantities  (in  thousands  of  acre-feet)  and  mean  EC 
values  (in  micromhos)  of  precipitation,  surface  water,  and  ground- 
water, plus  "additional  salt"  (in  tons)  applied  to  the  valley 
floor  drainage  area  of  each  reach.   The  total  inflowing  salt  load 


to  the  area  was  the  sum  of  the  products  of  the  above  flow  and  EC 
values,  expressed  in  tonsi/,  plus  the  tons  of  additional  salt 
added  by  fertilizers,  soil  conditioners,  etc. 

The  total  of  the  above  inflow  quantities  was  then  reduced  by  a 
depletion  quantity  which  consisted  of  the  consumptive  use  of  irri- 
gated crops  and  native  vegetation  plus  soil  moisture  requirements. 
The  result  was  the  net  accretion  of  flow  to  the  river  for  each 
month  (Figure  3)'   Assuming  all  inflowing  salts  were  returned  to 
the  stream,  the  EC  of  net  accretions  was  obtained  by  dividing  the 
total  inflowing  salt  load  as  computed  above  by  the  net  accretion 
flow.   By  inspection,  these  EC  values  were  then  modified  where 
they  deviated  significantly  from  known  historic  EC  values  of  irri- 
gation return  flow  or  runoff.   Final  flow  and  EC  results  became 
input  data  for  the  "other  accretions"  model  input  source  listed  on 
the  key  to  Figure  2.  Additional  study  is  required  on  the  accre- 
tion portion  of  the  model  so  that  it  will  more  consistently  produce 
more  realistic  results  of  both  flow  and  EC.   The  accretion  sub- 
routine is  explained  in  more  detail  in  "Computational  Methods "# 


1/  0.735  parts  per  million  (ppm)  of  total  dissolved  solids  (TDS) 
=  1  ton  of  TDS  per  1,000  acre-feet  of  water.   In  the  Sacramento 
River  0.735  ppm  is  roughly  equivalent  to  1  micromho  of  EC.  There- 
fore, 1  micromho  was  assumed  equal  to  1  ton  of  TDS  per  1,000  acre- 
feet  of  water. 
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Key  to  Figure  2,  "Water  Quality  Model  Input  Sources" 

Number 

Description 

Number 

Description 

1 

Shasta  Reservoir  Release 

47 

Additional  Losses 

2 

Spring  Creek  Tunnel 

50 

R.  D.  787(b)  Drain 

4 

Clear  Creek 

51 

Colusa  Basin  Drain 

5 

Cow  Creek 

55 

Sacramento  Slough 

6 

Cottonwood  Creek 

-  56 

Natomas  Cross  Canal 

7 

Battle  Creek 

■.■  57 

Additional  Inflow 

8 

Paynes  Creek   ,     , "  ' 

58 

M  &  I  Wastes 

9 

Redbank  Creek 

59 

Other  Accretions* 

10 

Additional  Inflow 

60 

Diversions 

11 

Other  Accretions* 

61 

Weir  Spills 

12 

Diversions 

62 

Additional  Losses 

13 

Additional  Losses 

66 

R.  D.  1000  Drain 

16 

Antelope  Creek 

67 

Natomas  Main  Canal 

17 

Elder  Creek 

68 

Natomas  E.  Main  Drain 

18 

Mill  Creek 

-  69 

Cache  Creek  (Eel  River 

19 

Thomes  Creek 

Diversion,  Southerly 

20 

Deer  Creek 

Route) 

21 

Additional  Inflow 

70 

Additional  Inflow 

22 

Other  Accretions* 

71 

M  &  I  Wastes 

23 

Diversions 

72 

Other  Accretions* 

24 

Weir  Spills 

73 

Diversions 

25 

Additional  Losses 

74 

Weir  Spills 

28 

Stony  Creek  (Eel  River 

75 

Additional  Losses 

Diversion,  Northerly 

79 

Additional  Inflow 

Route) 

80 

M  &  I  Wastes 

29 

Big  Chico  Creek 

81 

Other  Accretions* 

30 

Additional  Inflow 

82 

Diversions 

31 

M  &  I  Wastes 

83 

Additional  Losses 

32 

Other  Accretions* 

86 

Oroville  Reservoir  Release 

33 

Diversions 

87 

Yuba  River 

34 

Weir  Spills 

88 

Bear  River 

35 

Additional  Losses 

89 

Other  Accretions* 

36 

Butte  Creek 

90 

Diversions 

11 

Butte  Slough 

91 

Additional  Losses 

R.  D.  70  Drain 

94 

Folsom  Reservoir  Release 

41 

R.  D.  108  Drain 

95 

M  &  I  Wastes 

42 

R.  D.  787(a)  Drain 

96 

Other  Accretions* 

i2 

Additional  Inflow 

11 

Diversions 

Other  Accretions* 

Additional  Losses 

45 

Diversions 

46 

Weir  Spills 

*  Comprised  of  final  flow  and  EC  results. 
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COMPUTATIONAL  METHODS 


The  mean  monthly  EC  at  each  river  station  was  assumed  equal  to 
the  flow-weighted  mean  monthly  EC  of  all  Inflows  and  outflows  In 
the  reach  under  consideration.   The  products  of  each  flow  and  its 
EC  were  summed  algebraically  and  divided  finally  by  the  total 
flow  at  the  downstream  station.  This  may  be  expressed  as: 

Cy  =  (QxCx  +  QtiCti  + Qm^Tn  +  QaCQ  +  Vm  +  ^a^a  - 

QjjCjj   -  QgCs    -   Ql^l)  /   (%  +     QtI   + Qrn  +  Qq  +  ^M  + 

^A   -   ^D   -  %    -  <^l) 

Eq.    (1) 

where, 

C  =  mean  monthly  EC  in  micromhos 

Q  =  monthly  flow  in  1,000  acre-feet 

Y  =  downstream  station  In  reach 

X  =  upstream  station  in  reach 

T  =  measured  tributary 

n  =  number  of  tributaries 

G  «  estimated  groundwater  or  other  inflow  not  accounted 

for  in  the  other  model  inputs 
M  =  municipal  and  Industrial  waste  discharges 
A  =  unmeasured  accretions  (irrigation  return  flow  and 

additional  surface  water  runoff) 
D  =  diversions  from  the  river  for  beneficial  use 
S  =  spillage  from  the  river  over  flood  control  weirs 
L  =  unmeasured  losses  from  the  river 

The  numerator  in  Eq.  (l)  is  equivalent  to  the  total  salt  load 
fsSL)  in  tons  entering  the  river  in  the  reach  under  consideration 
(the  product  of  flow  in  1,000  acre-feet  and  EC  in  micromhos  is 
approximately  equal  to  tons  of  salt).   The  denominator  in  Eq.  (l) 
is  the  flow  (SQ)  in  the  river  at  the  downstream  station  in 
1,000  acre-feet.   Eq.  (l)  then  may  be  expressed  more  simply  as 
Cy  =  ESL/SQ. 

Plow  quantities  and  EC  values  for  item  A,  unmeasured  accretions, 
in  Eq.  (l)  were  computed  by  the  accretion  subroutine.   As  an 
option,  however.  Individual  accretion  values  could  be  specified 
in  the  input  data,  in  which  case  the  computed  values  would  not 
be  used.   The  accretion  subroutine  is  illustrated  schematically 
in  Figure  3  and  mathematically  below. 
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Solving  for  accretion  flow; 

Qa  =  R(Qp  +  Qs  +  Qg  -  "iET)  Eq.  (2) 

where, 

Q^  =  net  accretion  flow  In  1,000  acre-feet 

R  =  coefficient  of  runoff 

Qp  =  adjusted  precipitation  in  1,000  acre-feet 

Qs  =  applied  surface  water  In  1,000  acre-feet 

Qq  =  applied  ground  water  In   1,000  acre-feet 

Qgrp  =  evapotransplratlon  3n  1,000  acre-feet 

Solving  for  accretion  EC; 

Ca  =  P(QpCp  +  QsCs  +  QqCg  +   As)/Qa  Eq.    (3) 

where, 

Ca  =  EC  of  net  accretion  In  mlcromhos 
F  =  salt  return  factor 

Cp*  ^s,  Cq  =  EC  in  mlcromhos  of  precipitation,  applied 
surface  water  and  applied  groundwater, 
respectively 
As  =  additional  salt  in  tons  (fertilizers,  soil  condi- 
tioners, etc.) 

The  coefficient  of  runoff,  R,  has  not  yet  been  established  with 
any  reliability  for  subareas  of  the  Sacramento  Valley.   It  was  not 
used  in  the  manner  indicated  by  Bq.  (2)  for  the  1969-7O  study. 
Instead,  storm  water  runoff  and  irrigation  return  flow  values  were 
modified  by  inspection  where  they  deviated  significantly  from 
measured  historic  values. 

Coefficients  should  be  developed,  however,  so  that  the  model  can 
more  rationally  simulate  the  processes  which  produce  storm  water 
runoff  and  return  flow  from  applied  irrigation  water. 

The  salt  return  factor,  P,  as  in  the  case  of  the  flow  factor,  R, 
requires  additional  study  for  Its  development.   It  is  directly 
related  to  quantities  of  applied  salts  retained  in  the  soil 
during  the  irrigation  season  and  to  quantities  of  salts  leached 
from  the  soil  during  periods  of  heavy  rainfall. 

Base  Period  for  Hydrologlc  Applications 

The  hydrologlc  base  for  the  I965  model  was  the  period  from  1953 
to  1961,  which  covered  a  wide  range  of  hydrologlc  conditions. 
For  development  of  the  current  model,  time  and  budgetary  limita- 
tions restricted  selection  for  data  processing  to  only  three 
hydrologlc  years--1965,  1961,  and  1955. 

The  1965  water  year  was  chosen  for  preliminary  development  of  the 
model's  structure  and  for  writing  the  computer  program  because  it 
represented  an  approximate  present  level  of  development.   Also, 
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an  expanded  source  of  water  quality  data  was  readily  available 
from  the  Department's  data  collection  programs. 

For  1965  conditions,  the  objective  was  to  devise  a  model  which 
would  reasonably  simulate  the  Sacramento  Valley  water  resource 
system.   The  main  concern  was  to  develop  a  method  for  computing 
flow  and  EC  of  valley  floor  accretions. 

Net  accretions  (and  depletions)  first  were  approximated  by  evalu- 
ating differences  in  flow  and  EC  between  the  upstream  and  down- 
stream stations  in  each  reach,  based  on  recorded  values  for  the 
river  stations.   That  method  was  not  satisfactory,  however, 
because  net  accretions  did  not  indicate  individual  accretion 
sources,  such  as  irrigation  return  flow,  storm  water  runoff, 
groundwater  Inflow,  etc.   The  nature  of  certain  losses  in  flow, 
other  than  recorded  diversions  and  weir  spills,  also  could  not  be 
accounted  for.   A  more  refined  accretion  analysis  method  was 
required.   The  accretion  subroutine  was  the  result.   This  was 
used  in  subsequent  runs  of  the  model,  with  the  1961  water  year  as 
the  base  hydrology  for  further  refinements. 

The  1961  water  year  was  selected  because  the  Sacramento  River 
Water  Pollution  Survey  conducted  at  that  time  provided  an  addi- 
tional source  of  data,  including  data  on  agricultural  drains  not 
before  or  since  Included  in  the  Department's  data  collection 
programs . 

To  develop  a  better  rationale  for  use  in  projections  of  flow  and 
EC  under  future  levels  of  development,  flow  and  EC  values  for 
1961  were  first  computed  for  each  river  station,  using  the  accre- 
tion subroutine  for  accretions,  and  results  were  then  compared 
with  recorded  values.   Where  deviations  between  computed  and 
recorded  values  were  excessive,  adjustments  were  made  in  either 
valley  floor  inflow  or  losses  (provided  the  adjustments  could  be 
Justified  on  some  rational  basis  or  logical  assumption). 

The  1953  water  year  was  used  to  verify  the  model  because  it  was 
the  second  year  for  which  comprehensive  diversion  data  were 
available  from  the  Trial  Distribution  studies  of  flow  in  the 
Sacramento  Valley  system  conducted  by  the  U.  S.  Bureau  of  Recla- 
mation, the  Department  of  Water  Resources,  and  local  water  users. 
Water  quality  data  also  were  available  from  the  Department's 
data  programs,  which  began  systematic  collection  of  surface  and 
groundwater  quality  data  in  1951.  Results  of  the  1955  verifica- 
tion run  are  discussed  under  "Model  Verification". 
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6 J   Net  accretions  (including  natural  salt  loadings) 

7)  Municipal  and  industrial  waste  discharges 

8)  Additional  inflow  to  the  river 

9)  Direct  diversions  from  the  river  for  beneficial  use 
Spillage  from  the  river  over  flood  control  weirs 


11)  Additional  losses  of  river  flow 

Tables  1  and  2  are  copies  of  machine  printouts  which  show  how  the 
above  factors  were  incorporated  into  the  model.   Table  1,  "Analysis 
of  Accretions  for  Reach  from  Below  Colusa  Basin  Drain  to  Verona", 
!  shows  those  factors  which  were  considered  in  the  model's  accretion 
I  subroutine  and  which  contributed  to  valley  floor  accretions  of 
j  both  flows  and  salts  in  each  reach  of  the  river.   The  particular 
I  reach  considered  in  Table  1  is  Reach  G,  below  Colusa  Basin  drain 
to  Verona.   Net  accretions  contain  agricultural  return  flows,  run- 
off from  precipitation,  or  a  combination  of  the  two,  depending 
upon  the  particular  month  under  consideration. 

Table  2,  "Estimated  I96I  Monthly  Flow  and  Mean  Conductivity  in 
Sacramento  River  System,  Below  Colusa  Basin  Drain  to  Verona", 
shows  all  of  the  inflows  (including  previously  computed  accre- 
tions) and  outflows  which  occur  within  the  same  reach.   The  final 
values  are  net  flows  and  conductivities  at  Verona,  the  downstream 
station  of  the  reach.   Zero  values  are  shown  for  Sacramento  Slough 
and  Natomas  Cross  Canal  because  these  are  channels  for  agricul- 
tural drainage  or  other  surface  water  runoff,  values  for  which,  in 
this  case,  have  been  computed  in  the  accretion  analysis  and  are 
included  in  the  item  "other  accretions".   If  desired,  flow  and  EC 
values  may  be  specified  for  these  two  inputs  rather  than  computing 
them  in  the  accretion  analysis.   In  that  manner,  the  effects  on 
river  EC  of  any  individual  agricultural  drain  at  any  point  in  the 
system  can  be  evaluated  by  the  model.   The  option  of  specifying 
accretion  values  is  discussed  further  under  "Net  Accretions". 

Precipitation 

Historic  monthly  quantities  of  precipitation  falling  on  the  valley 
floor  area  of  each  reach  were  available  through  the  Department's 
studies  (23)  of  consumptive  use  and  water  requirements.   In  those 
studies,  estimated  average  monthly  rainfall  was  computed  by  multi- 
plying measured  rainfall  at  key  stations  by  an  appropriate  Index 
factor  to  yield  the  rainfall  on  the  given  land  area  in  terms  of 
thousands  of  acre-feet. 

Precipitation  quantities,  before  being  used  as  inputs  to  the  water 
quality  model,  were  modified  by  the  effects  of  soil  moisture  stor- 
age and  deep  percolation  to  groundwater.  Durlrig  some  months,  a 
portion  of  direct  precipitation  served  to  replenish  soil  moisture 
storage  and/or  storage  in  the  groundwater  basin,  and  precipitation 
values  were  adjusted  by  deducting  those  quantities.   During  other 
months,  water  was  withdrawn  from  soil  moisture,  and  precipitation 
values  were  adjusted  by  adding  the  amounts  withdrawn.   The  latter 
condition  was  obtained  generally  during  the  irrigation  season  and 
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TABLt  1.   ANALYSIS   or   ACCKCTIOnS 
TOR  RMOH  PROM 
BCLO»    ChO    to    VCRONA 


1961  WATCT  TWR 


ocrnaEf)  NOvExeE*)  ueccmbcr   jANUAxr  revxuARr    march 


AUSUtT    St^TCMSCa 


•••  eONTRlHUTIONJ  TO  ACCHETION  ••• 


PHECIPITATIDN 


FLO*         15 
t.C.        100 

SALT        1500 

116 

100 
11600 

1«0 
4600 

«0 

100 
9000 

9S 

104 
9500 

I3T 
100 
liTOO 

9« 

lOS 
9400 

69 

100 
69«0 

20 

100 
2000 

i 

101 
500 

t 
100 
»00 

13 

loo 

1300 

SURFACE  WATtN 
FLU*          6 
t.C.        1B9 
SALT        1135 

1 
IVb 
1"»5 

1 
IHO 

180 

0 
100 

0 

0 
157 

0 

0 
16T 

0 

61 
207 
12611 

116 

26M 
31036 

112 
259 
2B9S« 

120 
2»' 
29665 

111 

221 

24^6* 

•  40 

197 
9471 

GHOUNR  KATEH 

FLO*         ?<> 
t.C.       ton 
SALT       ?6100 

0 
900 

0 

0 
»00 

0 

0 
900 

0 

0 

900 

0 

0 
900 

0 

13 

900 
11700 

10 
900 
9000 

20 

900 

tiooo 

19 

900 
17100 

20 
900 
ItlOOO 

* 

900 
7200 

AOOITIOMAL  SALT 
TOTAL  SALT 

?a7 

OtPLETlON  FLOW 


NET  ACCBFTIUNS 
FLOW  0 

E.C.  0 

SALT 


4780 

9000 

9500      12706      33721 

46942 

48961 

34 

31 

66         95        16a 
•••  RESULTANT  ACCRETION  ••• 

178 

137 

U 
6a 

4760 

59 
350* 
20650 

29        32         0 

32a        397          0 
9500      12706          0 

17 
400* 
6800 

15 

400» 
6000 

43067 
126 


(•)   INOICATES   U   OR   C   VALUE  SPECIFIED  IN  INPUT  OATA.  COHPUTEO  VALUE  NOT  USED. 


TABLE  2.   ESTIMATED  I96I  MONTHLY  FLOW*  AND  MSAN  COKDUCTIVITY*  IN  SACRAMENTO  RIVER  SYSTEM 
RELOH  CBO  TO  VEMONA 
1961  MYOHOLOOT 


THIHUIArilFS  AND  OTmER 

INFLO.S 

OCT 

NOV 

DEC 

JAN 

'ER 

MAR 

APR 

MAT 

JUN 

JUL 

AUG 

SEP 

oELo*  cnn 

l-LO* 

361 

447 

731 

467 

1'45 

1316 

546 

429 

413 

426 

484 

4»1 

CONOUCTIVIIr 

ia9 

195 

lao 

lao 

1»7 

167 

207 

268 

?5V 

247 

221 

1V7 

FEATHFH-NlCOLAUa 

fLO* 

57 

157 

264 

IHO 

SOB 

489 

404 

281 

108 

40 

40 

32 

ClINIlUCtlVlTy 

147 

144 

120 

126 

105 

96 

90 

«9 

105 

133 

132 

IS' 

SACRAMtnro  SL 

ao* 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CONDUCTIVITY 

U 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

NATOMAS  A  CANAL 

►  LO* 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CONDUCTIVITY 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

AODL  INFLOW 

►  lo* 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

45 

CONDUCTIVIIY 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

n 

900 

M  ANO  1  WAiteS 

►  LOW 

0 

u 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CONDUCTIVITY 

700 

700 

700 

700 

700 

700 

700 

700 

700 

700 

700 

700 

OTMER  ACCRETIONS 

Flo* 

0 

55 

13 

59 

?9 

32 

0 

17 

15 

0 

10 

0 

CONOUCTIVITY 

0 

214 

36S 

350 

328 

397 

0 

400 

400 

0 

400 

0 

DIVERSIONS  ANO  OT"£R 

OUTKLOlS 

DIVERSIONS 

►  LO* 

0 

0 

0 

0 

IT 

IS 

18 

19 

CONOOCnviIY 

0 

0 

0 

0 

207 

268 

259 

247 

221 

197 

*EIR  SPILLS 

FLO* 

0 

0 

0 

0 

0 

0 

0 

0 

CONDUCTIVITY 

1B3 

184 

166 

180 

146 

153 

157 

200 

231 

237 

2U 

257 

OTHER  LUiStS 

FLO* 

0 

0 

0 

0 

0 

0 

0 

0 

CONDUCTIVITY 

0 

0 

0 

0 

0 

0 

0 

0 

NET  VALUES  AT 

VERONA 

►  LO* 

41» 

659 

1008 

T06 

1RI2 

IB39 

943 

710 

921 

448 

515 

512 

CONDUCT  1 VI  FY 

1»3 

184 

166 

180 

146 

153 

157 

200 

231 

237 

2ia 

257 

•  Monthly  flow  1»  ihoan  In  1,000  acre-feet  end  mean  conductivity  ia  ahovn  In  nlcronhoa. 
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often  would  indicate  that  direct  precipitation  had  occurred  when 
in  fact  it  was  a  consumptive  use  of  stored  precipitation  that  had 
occurred  earlier. 

No  systematic  records  were  available  on  the  mineral  content  of 
rainfall  in  the  Sacramento  Valley.   Chemical  analyses  of  snow  on 
the  west  slope  of  the  Sierra  Nevada  (13)  have  shown  total  dis- 
solved solids  ranging  from  less  than  1  to  11  parts  per  million, 
with  a  mean  value  of  about  5  ppm.   In  this  study,  EC  of  rainfall 
was  initially  assumed  to  be  a  constant  I5  micromhos,  corresponding 
to  about  10  ppm  of  total  dissolved  solids.   In  later  model  runs, 
the  EC  value  was  increased  from  I5  to  100  micromhos  to  reflect  the 
mineral  pickup  of  rainwater  as  it  flowed  over  the  soil  and  rocks 
of  the  valley  floor.   This  produced  a  more  realistic  result,  since 
it  corresponded  more  closely  to  that  of  surface  water  runoff  into 
the  river. 

Applied  Surface  Water 

The  term  "applied  surface  water"  as  used  in  this  study  refers  to 
amounts  of  water  diverted  from  the  Sacramento  River  or  its  tribu- 
taries and  used  to  irrigate  crops  in  the  area  under  consideration. 
These  are  not  necessarily  the  same  amounts  of  water  per  reach  as 
those  included  in  "diversions",  which  will  be  discussed  later. 
Water  diverted  from  the  river  in  one  reach  may  be  applied  for  irri- 
gation in  another  area  which  contributes  return  flow  to  a  reach 
downstream.   For  example,  the  Glenn-Colusa  Canal  carries  water 
diverted  from  Reach  "B"  and  applied  in  Colusa  Basin  where  drain- 
age returns  to  Reach  "F". 

For  any  historic  year,  monthly  diversion  flow  quantities  are 
recorded  in  publications  on  surface  water  flow  (16,  19).   These 
recorded  values,  where  applicable,  were  used  as  the  basic  sur- 
face water  supply  in  accretion  analyses  for  each  reach.   In  some 
areas  unmeasured  amounts  of  surface  water  also  are  diverted  from 
minor  tributaries  crossing  the  area.   In  such  instances  measured 
diversion  quantities  alone  may  be  insufficient  for  use  as  the 
total  surface  water  supply.   Where  estimated  depletions  in  a 
particular  area  exceeded  the  estimated  total  water  supply  (precip- 
itation, measured  stream  diversions,  and  estimated  groundwater 
applications)  the  deficit  was  assumed  to  be  made  up  by  unmeasured 
diversions  from  minor  tributaries. 

The  EC  of  applied  surface  water,  in  general,  was  assumed  to  be 
the  same  as  the  EC  of  the  Sacramento  River  at  the  upstream  sta- 
tion of  the  reach  from  which  the  water  was  diverted.   For  Colusa 
Basin,  however,  the  average  of  EC's  at  Colusa  and  Hamilton  City 
was  used. 

Applied  Groundwater 

Very  little  data  are  available  on  measured  quantities  of  ground- 
water applied  historically  in  the  Sacramento  Valley.   Preliminary 


19 


estimates  for  use  in  the  model  were  provided  by  the  Coordinated 
Statewide  Planning  Program  of  the  Northern  District  office  of  the 
Department. 

EC's  of  groundwater  used  in  the  model  were  the  approximate  areal 
averages  of  EC  data  gathered  from  wells  since  1951  "by  the  Depart- 
ment's groundwater  quality  data  program  (l6,  17). 

Additional  Salts 

The  total  quantities  of  salts  in  a  given  area,  in  addition  to  those 
occurring  there  naturally  and  those  brought  in  by  inflowing  water, 
include  additional  salts",  those  added  by  fertilizers,  soil  con- 
ditioners, etc.   Quantities  of  these  added  salts  for  use  as  model 
inputs  were  obtained  from  various  sources.   Available  publications 
by  the  University  of  California  Experiment  Stations  were  reviewed 
and  additional  information  was  obtained  by  interviewing  farm 
advisors  and  other  knowledgeable  individuals. 

The  amounts  of  additional  salts  are  minor  in  relation  to  those 
carried  by  the  river  and  its  tributary  inflows.   For  example,  in 
1961  a  total  of  59  tons  of  salts  were  estimated  to  have  been 
returned  to  the  river  from  such  sources  as  fertilizers,  soil  amend- 
ments, etc.,  during  April,  the  month  of  greatest  application.   The 
total  salt  loading  in  the  river  at  Freeport  during  that  month  was 
estimated  at  163,300  tons.   These  additional  salts  were  less  than 
0.04  percent  of  the  total  salts  in  the  river.   In  terms  of  nutri- 
ents, however,  fertilizers  may  be  of  much  greater  significance. 

Depletions 

Monthly  depletions  in  water  supplies  to  meet  the  consumptive  de- 
mands of  agricultural,  urban  or  native  lands  could  not  be  computed 
directly  by  the  water  quality  model.  A  side  model,  programmed  for 
digital  computer,  was  used  to  make  these  computations.   This  side 
model  was  developed  by  the  Department  as  part  of  the  continuing 
Joint  DWR-USBR  Central  Valley  operation  studies  (23),  which  de- 
veloped a  program  commonly  referred  to  as  the  C.U,-2  Program. 

The  C.U.-2  Program,  which  will  simulate  actual  field  moisture  con- 
ditions, is  capable  of  handling  the  lag  in  time  between  the  time 
rainfall  or  irrigation  water  enters  the  soil  reservoir  and  the 
time  it  is  lost  to  evapotranspiration.   Water  stored  in  the  soil 
is  treated  as  part  of  the  water  supply  in  the  month  in  which  it 
is  used.   Depletions  in  any  given  month  are  essentially  restricted 
to  evapotranspiration. 

The  C.U.-2  Program,  simulating  the  action  of  the  soil  reservoir, 
stores  precipitation  that  falls  in  the  winter  months  for  use  in 
the  spring  and  summer.  Applied  irrigation  water  also  is  retained 
in  the  soil  or  released  in  the  same  manner,  but  no  water  is 
stored  in  the  winter.   With  the  C.U.-2  Program  it  is  possible  to 
consider  the  differences  in  rooting  depth  and  growing  seasons  of 
various  crops,  the  moisture  retention  of  various  soil  textures, 
or  the  addition  of  ponding  water  for  rice  during  April  and  May. 
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For  this  study,  the  C.U.-2  model  was  used  to  compute  both  monthly 
and  annual  depletions  for  eight  major  agricultural  crops  and  for 
urban  and  native  lands.   These  land  uses  and  acreages  for  the 
selected  water  years  are  shown  in  Appendix  C. 

Net  Accretions 

Accretions  such  as  storm  water  runoff  and  return  flows  from  irri- 
gation are  computed  in  the  accretion  subroutine  of  the  model, 
unless  specified  directly  in  the  input  data.   Examples  of  flow 
and  EC  values  of  computed,  as  well  as  specified  net,  accretions 
are  shown  on  Table  1,  with  specified  values  designated  by  an  aster- 
isk.  Final  net  accretion  values  are  carried  as  other  accretions" 
for  each  reach,  as  shown  on  Table  2. 

In  the  accretion  subroutine,  all  salts  applied  to  an  area  during 
a  particular  month  were  assumed  returned  to  the  stream  that  same 
month.   This  is  not  necessarily  true;  and  strong  evidence  exists 
to  show  that,  in  the  Sacramento  Valley,  some  of  those  salts  gen- 
erally are  retained  in  the  area  during  the  irrigation  season  and 
are  returned  to  the  stream  when  leached  from  the  soil  during 
periods  of  heavy  rainfall.   For  that  reason,  the  option  of  speci- 
fying the  flow  and/or  EC  of  accretions  was  incorporated  into  the 
structure  of  the  model.   Whenever  an  accretion  value  is  specified, 
the  computed  value  will  not  be  used  as  input  to  the  model.   Observ- 
able patterns  of  historic  accretions  can  thereby  be  reproduced. 
During  the  brief  period  allotted  to  this  study,  a  detailed  analy- 
sis of  these  possible  patterns  could  not  be  undertaken.   However, 
available  flow  and  EC  data  on  return  flows  and  runoff  were  in- 
dispensable guides  to  probable  magnitudes  of  values  specified. 

Accretion  flow  values  were  specified  most  frequently  for  months 
of  relatively  heavy  rainfall.   For  those  periods,  the  accretion 
subroutine  often  produced  net  accretion  flow  values  considerably 
in  excess  of  apparent  historical  quantities  of  runoff.   In  the 
model,  more  realistic  flow  results  were  obtained  when  runoff  in 
the  Sacramento  Valley  floor  during  any  one  month  did  not  exceed 
30  percent  of  that  month's  total  of  precipitation,  applied  sur- 
face water,  and  applied  groundwater.   For  that  reason,  net  accre- 
tion values  were  limited  to  30  percent  of  precipitation  and 
applied  water.   Admittedly,  this  30  percent  factor  was  a  some- 
what arbitrary  runoff  coefficient.  A  more  valid  coefficient 
undoubtedly  would  vary,  depending  upon  antecedent  precipitation 
and  soil  moisture  conditions,  types  of  terrain,  and  other  factors. 
Further  refinement  of  the  accretion  analysis  portion  of  the  model 
should  produce  more  valid  runoff  coefficients. 

Municipal  and  Industrial  Waste  Discharges 

The  analysis  of  accretions  (Table  1)  provided  flow  and  EC  esti- 
mates of  agricultural  drainage  and  surface  water  runoff,  but  did 
not  include  estimates  of  municipal  and  industrial  (M&I)  wastes 
discharged  to  the  river  system.   These  wastes  were  treated  as  an 
individual  model  input,  as  shown  on  Table  2. 
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For  the  1955-62  period  covered  by  DWR  Bulletin  No.  68-62  (6), 
quantities  of  M&I  wastes  discharged  directly  to  the  Sacramento 
River  and  its  tributaries  varied  from  62,330  to  79,230  acre-feet 
per  year.l/  Since,  in  general  the  bulletin  considered  only  indi- 
vidual waste  discharges  exceeding  500  acre-feet  per  year,  the 
above  total  quantities  could  be  low.  For  the  1955  verification 
run,  a  total  of  84,000  acre-feet  was  used,  with  60,000  acre-feet 
in  the  Sacramento-Freeport  reach;  12,000  acre-feet  between 
Sacramento  and  Colusa  (introduced  in  the  Verona -Sacramento  reach); 
and  12,000  acre-feet  above  Colusa  (introduced  in  the  Hamilton 
City-Colusa  reach). 

EC  of  M&I  wastes  during  the  1955-62  period  at  the  City  of  Sacramento 
sewage  treatment  plant  varied  from  642  to  738  micromhos,  with  an 
average  of  699  micromhos.  An  EC  of  700,  therefore,  was  considered 
representative  for  M&I  wastes  and  was  used  as  a  constant  throughout 
the  system  for  both  historic  and  future  levels  of  development. 
Sewage  treatment  processes  by  1990  in  the  Sacramento  Valley  were 
not  expected  to  change  significantly  with  respect  to  EC;  therefore, 
a  total  dissolved  solids  in  M&I  wastes  were  assumed  to  remain  about 
the  same  as  they  are  now.  However,  nutrients  and  toxic  materials 
contained  in  M&I  wastes  are  more  likely  to  adversely  affect  river 
water  quality  than  would  total  dissolved  solids. 

Additional  Inflow 

In  developing  the  current  model  an  essential  step  was  to  make  an 
approximate  hydrologic  balance  between  upstream  and  downstream 
stations  of  each  reach  for  which  recorded  flow  data  were  available. 
After  accounting  for  all  known  or  estimated  inflows  and  outflows, 
significant  additional  inflows  were  needed  to  achieve  flow  balance 
in  some  reaches  during  parts  of  both  the  196I  and  1965  water  years. 

A  typical  example  may  be  observed  in  Table  2.   September  I96I 
shows  that  an  additional  inflow"  of  45,000  acre-feet  (TAF)  was 
introduced  (arbitrarily)  in  the  reach  below  CBD  to  Verona.  The 
analysis  of  accretions  (Table  1)  shows  zero  net  accretions  for 
September,  which  appears  reasonable  in  view  of  the  relative 
amounts  of  precipitation,  applied  ground  and  surface  water,  and 
depletion  flows.   The  difference  between  recorded  and  computed 
flows  at  Verona  was  56  TAF  (523  TAF  versus  467  TAF,  respectively) 
before  adding  the  45  TAF  additional  inflow.  Although  the  adjusted 
value  at  Verona  was  still  less  than  the  recorded  value  by  11  TAF, 
the  new  computed  value  at  Sacramento  was  brought  to  within  4  TAF 
of  the  recorded  value  with  no  further  adjustments.   The  rela- 
tively small  difference  of  4  TAF  should  be  within  the  probable 
error  of  flow  measurements. 

After  the  probable  amount  of  additional  inflow  has  been  deter- 
mined, the  most  likely  source  of  such  inflow  must  be  considered 


1/  Additional  data  on  M&I  waste  discharge  are  shown  in 
15WR  Bulletin  No.  Ill,  bibliography  entry  No.  26. 
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so  that  an  appropriate  EC  value  can  be  assigned  to  it.   On  the  basis 
of  EC,  these  additional  inflows  could  most  often  be  best  attributed 
to  groundwater.   Assignment  of  known  groundwater  EC  values  to  the  in- 
flow (as  was  done  for  the  45  TAF  additional  inflow  shown  on  Table  2) 
and  operation  of  the  model  could  approximately  reproduce  recorded 
EC  values  at  the  river  stations. 

Where  hydrologic  conditions  are  favorable,  groundwater  inflow  un- 
doubtedly occurs  along  the  river.   When  the  river  stage  is  rising, 
water  will  be  recharged  into  bank  storage.   When  the  stage  drops, 
water  will  tend  to  be  withdrawn  from  bank  storage  to  the  river. 
Wherever  the  adjacent  groundwater  level,  whether  in  bank  storage  or 
in  the  groundwater  basin,  becomes  higher  than  the  surface  of  the 
river,  groundwater  inflow  can  occur.   No  firm  criteria  were  devel- 
oped in  this  study  upon  which  to  base  the  occurrence  or  quantities 
of  groundwater  inflow.   Historically,  the  hydrologic  balance  served 
as  a  guide. 

The  EC  of  groundwater  inflow  generally  was  assumed  to  be  the  same 
as  that  measured  in  the  adjacent  groundwater  basin.   Wherever  in- 
flow appeared  to  be  derived  mainly  from  bank  storage,  EC  was  as- 
sumed to  be  considerably  lower  than  the  EC  of  water  pumped  from 
wells,  particularly  in  those  areas  containing  more  highly  miner- 
alized groundwater. 

Groundwater  inflow  from  subsurface  agricultural  drainage  can  affect 
the  quality  of  Sacramento  River  water.   A  recent  study  (15)  reported 
by  the  University  of  California  Agricultural  Extension  Service  in 
Davis  was  an  important  guide  in  judging  the  probable  significance  of 
the  water  quality  problems  caused  by  such  inflow.   Its  findings  as- 
sisted in  assigning  flow  and  EC  values  to  model  inputs.   The  study 
estimated  the  extent  to  which  future  subsurface  agricultural  drain- 
age will  be  a  problem  in  Tehama,  Glenn,  Colusa,  and  Yolo  Counties  in 
20'^0.   Findings  for  the  four-county  area  are  summarized  here, 

1)  Area  that  will  need  artificial  subsurface  drainage: 
75,150  acres 

2)  Quantity  of  subsurface  drainage  water:  232.2  cfs 
(13.9  TAF  per  month) 

3)  Salinity  (EC)  range:  5OO  to  5,000  micromhos 

Almost  70  percent  of  the  above  drainage,  or  159^7  cfs,  would  enter 
the  Colusa  Basin  drain  at  an  EC  ranging  from  500  to  5,000  micromhos. 
The  remaining  30  percent,  or  72.5  cfs,  would  enter  the  Sacramento 
River  at  Stony  Creek  or  above,  with  an  EC  ranging  from  5OO  to 
2,000  micromhos.   The  total  of  this  drainage  is  small  in  relation 
to  other  valley  floor  accretions  to  the  river,  amounting  to  only 
about  6  percent  of  the  projected  agricultural  drainage  that  would 
enter  the  river  between  Keswick  and  Knights  Landing. 

Diversions 

Direct  diversions  from  the  river  in  a  particular  reach  may  not  be 
the  same  as  diversions  for  "applied  surface  water"  in  that  same 
reach,  because  a  portion  of  the  diversions  may  be  applied  in  the 
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area  of  a  different  reach.   For  historic  conditions,  data  on  diver- 
sions of  streamflow  are  available  from  state  publications  on  sur- 
face water  flow  (l6,  19).   For  future  conditions,  diversion 
quantities  were  related  to  quantities  of  applied  surface  water 
which  would  be  withdrawn  from  the  river,  but  not  necessarily  ap- 
plied in  a  particular  reach. 

The  EC  of  diversions  was  assumed  to  be  the  same  as  the  EC  in  the 
river  at  the  upstream  station  of  each  reach.   For  prediction  model 
runs,  the  upstream  EC  would  be  the  only  "known"  EC  river  value  in 
the  reach,  since  it  would  have  been  the  final  result  of  computa- 
tions for  the  reach  immediately  upstream. 

Weir  Spills 

Five  weirs  along  the  Sacramento  River  divert  floodflows  from  the 
river  into  bypass  channels.   These  are  the  Moulton,  Colusa,  Tisdale, 
Fremont,  and  Sacramento  Weirs.   Historic  flow  over  these  weirs  is 
recorded  in  reports  on  surface  water  flow  (l6,  19,  12).  As  in 
other  diversions,  EC  of  weir  spill  flow  is  the  same  as  that  of  the 
river  at  the  weir  location.   With  one  exception,  this  is  consid- 
ered in  the  model  to  be  the  EC  at  the  upstream  station  in  each 
river  reach  under  consideration. 

The  exception  for  determining  the  EC  of  weir  spillage  occurs  in 
Reach  "G  ,  below  Colusa  Basin  Drain  (Station  10?)  to  Verona 
(Station  IO9).   In  that  reach,  the  computed  EC  at  Station  10? 
could  be  a  considerably  inaccurate  basis  for  computing  the  EC  of 
outflow  over  Fremont  Weir.   That  is  because,  under  high  flow  con- 
ditions, Fremont  Weir  spillage  would  likely  include  large  Feather 
River  flows  at  EC  values  differing  widely  from  those  in  the 
Sacramento  River  upstream  from  the  weir.   The  EC  of  the  Fremont 
Weir  spill,  therefore,  was  computed  as  that  of  the  mixture  of  flows 
in  the  Sacramento  River  below  CBD  and  the  Feather  River  at  Nicolaus. 

Additional  Losses 

As  discussed  under  "Additional  Inflow",  a  hydrologic  balance  be- 
tween upstream  and  downstream  stations  was  attempted  for  each 
historic  year  considered.   During  some  months,  losses  in  river 
flow  were  detected  with  could  not  be  accounted  for  in  measured 
diversions.  Additional  losses  undoubtedly  do  occur  under  certain 
hydrologic  conditions.   When  the  river  stage  is  rising,  water  can 
be  expected  to  flow  from  the  river  into  bank  storage  and  possibly 
as  recharge  to  the  adjacent  groundwater  basin  or  seepage  onto 
adjacent  lands. 

No  firm  criteria  or  patterns  for  such  losses  were  developed  during 
the  course  of  this  study.   However,  observation  showed  that  the 
most  substantial  losses  usually  occurred  during  periods  of  sharp 
increases  in  river  flow.   In  the  model,  the  EC  of  these  losses  was 
considered  in  the  same  manner  as  was  the  EC  of  diversions:   equiv- 
alent to  EC  of  the  river  at  the  appropriate  upstream  locations. 
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MODEL  VERIFICATION 


To  verify  the  model,  computed  EC  values  were  compared  graphically 
and  numerically  with  prototype  values  at  selected  river  stations. 
Because  data  for  comparisons  were  limited  to  only  one  year  of 
record  and  because  of  uncertainties  in  the  record  itself,  an  in- 
volved statistical  approach  toward  verification  was  not  considered 
appropriate.   Verification,  therefore,  was  largely  a  matter  of 
individual  judgment  guided  by  knowledge  of  the  Sacramento  River 
system,  the  structure  of  the  model,  and  the  adequacy  of  data  used. 

As  discussed  earlier,  the  I965  water  year  was  used  for  developing 
and  programming  the  structure  of  the  model;  the  1961  water  year, 
for  refinements  and  preliminary  verification;  and  the  1955  water 
year,  for  final  verification.   Several  machine  runs  were  made  for 
1961  until  reasonable  agreement  was  achieved  between  computed  and 
measured  values.   Successive  adjustments  in  flow  and/or  EC  of  cer- 
tain unmeasured  input  data  were  necessary,  primarily  for  additional 
inflow  or  additional  losses. 

Certain  "boundary"  conditions,  of  course,  were  considered  fixed. 
These  included  all  measured  or  estimated  flow  rates  and  EC  values 
of  reservoir  releases,  tributary  inflows,  diversions,  precipita- 
tion, applied  irrigation  water,  consumptive  use,  and  municipal  and 
industrial  waste  discharges.   Certain  other  constraints  were  intro- 
duced, such  as  a  maximum  runoff  coefficient  of  30  percent  and  ad- 
herence to  approximate  historic  EC  values  for  irrigation  return 
flow.   The  latter  constraint  was  necessary  because  the  accretion 
subroutine  often  produced  results  which  deviated  greatly  from  EC 
values  observed  in  irrigation  drains. 

The  same  basic  assumptions  and  rationale  employed  for  I96I  and 
1965  were  used  for  the  1955  verification  run  and  computed  results 
were  compared  with  1955  recorded  data.   No  adjustments  were  made 
on  the  basis  of  recorded  EC  data  for  the  river  stations,  so  that 
the  model's  probable  accuracy  for  use  in  future  predictions  could 
be  evaluated. 

The  model  was  balanced  hydrologically ,  however,  against  recorded 
streamflow  rates  at  terminal  river  stations  in  each  reach  to 
account  for  unmeasured  accretions  or  depletions,  such  as  ground- 
water inflow  and  seepage  losses. 

Figures  4  through  8  show  comparisons  between  the  1955  recorded 
and  computed  EC  values  and  recorded  and  computed  river  flow  at  sta- 
tions where  historid  data  a  *e  available.  ALttiough  most  recorded  and 
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computed  flow  values  coincided,  they  were  not  necessarily  identical, 
but  deviations  were  not  great  enough  to  be  perceptible  on  the 
graphs.   The  close  correlation  in  flows  was  achieved  intentionally 
by  the  hydrologic  balance.  Additional  inflows  or  losses  in  each 
reach  were  adjusted  so  that  the  computed  flows  at  the  river  sta- 
tions would  approximate  the  measured  values.   For  future  conditions, 
the  model  was  balanced  hydrologically  by  reservoir  release  and 
river  flow  values  established  in  Central  Valley  water  project  op- 
eration plans. 
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Figure  4.   RECORDED  AND  COMPUTED  PLOW  AND  CONDUCTIVITY, 
SACRAMENTO  RIVER  AT  HAMILTON  CITY,  1955  WATER  YEAR 
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Figure   5.      RECORDED   AND   COMPUTED  FLOW   AND   CONDUCTIVITY, 
SACRAMENTO  RIVER  BELOW   COLUSA  BASIN  DRAIN    (KNIGHTS   LANDING), 
1955  WATER  YEAR 
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Figure  6.   RECORDED  AND  COMPUTED  PLOW  AND  CONDUCTIVITY,  FEATHER 
RIVER  AT  NICOLAUS,  1955  WATER  YEAR 
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Fleure   5.      RECORDED   AND   COMPUT^ED  FtOW   AND   CONDUCTIVITY, 
SACRAMENTO  RIVER  BELOW   COLUSA  BASIN  DRAIN    (KNIGHTS   LANDING), 
1955  WATER   YEAR 
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Figure  6.   RECORDED  AND  COMPUTED  PLOW  AND  CONDUCTIVITY,  FEATHER 
RIVER  AT  NICOLAUS,  1955  WATER  YEAR 
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Figure   7.      RECORDED   AND   COMPUTED  FLOW   AND  CONDUCTIVITY, 
AMERICAN   RIVER    AT   SACRAMENTO,    1955   WATER  YEAR 
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Figure   8.      RECORDED   AND   COMPUTED  FLOW   AND  CONDUCTIVITY, 
SACRAMENTO  RIVER   AT  SACRAMENTO,    1955  WATER  YEAR 
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Tables  3  through  5  list  EC  values  as  recorded  and  computed  at 
three  sites:   in  the  Sacramento  River  at  Hamilton  City,  in  the 
Sacramento  River  below  Colusa  Basin  drain,  and  in  the  Sacramento 
River  at  Sacramento.   Also  shown  are  deviations  of  computed 
values  from  recorded  values  for  each  month. 

At  Hamilton  -City  (Table  3)^  the  monthly  average  error  was  13«5  per- 
cent.  The  recorded  EC  value  for  December,  however,  definitely  did 
not  represent  the  monthly  mean, since  the  rate  of  flow  when  the 
sample  was  taken  was  two  times  as  great  as  mean  monthly  flow  in 
cubic  feet  per  second  (cfs).   Omitting  December,  the  monthly  aver- 
age error  was  7.5  percent. 


TABLE  3.   ERROR  ANALYSIS,  1955  WATER  YEAR 
Sacramento  River  at  Hamilton  City 


Conductivity 

in  Micromhos 

Deviation 

Month 

Recorded!./' 

Computed 

Micromhos 

Percent 

OCT 

121 

118 

-  3 

2.5 

NOV 

97- 

159 

+24 

17.8 

DEC 

174 

+77 

79.^ 

JAN 

130 

157 

+27 

20.8 

FEB 

158 

151 

-  7 

4.4 

MAR 

161 

138   ■   ■■ 

-23 

14.3   " 

APR 

144 

135 

-  9 

6.3 

MAY 

134 

128 

-  6 

4.5 

JUN 

133 

128 

-  5 

3.8 

JUL 

129 

127 

-  2 

1.6 

AUG 

133 

125 

-  8 

6.0 

SEP 

127 

128 
Total 

+  1 

0.8 

+66 

Avera 

ee  per  month 

+  5. 

5 

13.5 

Omitting  Dece 

mber:   Total 

-11 

Revised  Avera 

ge  per  month 

-  0. 

9 

7.5 

1/  Recorded  values  are  results  of  grab  samples  (composites  of  2) 
taken  on  only  one  day  each  month,  and  do  not  necessarily  repre- 
sent monthly  mean  values. 

2/  Sample  for  December  taken  during  peak  flow  (2  x  monthly  mean  ) 
and  EC  does  not  represent  monthly  mean. 
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TABLE  k,      ERROR  ANALYSIS,  1955  WATER  YEAR 
Sacramento  River  Below  Colusa  Basin  Drain 


Conductivity- 

in  Micromhos 

Deviation 

Month 

Recorded 

Computed 

Micromhos 

Percent 

OCT 

195 

151 

-44 

22.6 

NOV 

174 

188 

+14 

8.0 

DEC 

169 

227 

+58 

34.3 

JAN 

189 

190 

+  1 

0.5 

b'EB 

213 

181 

-32 

15.0 

MAR 

207 

180 

-27 

13.0 

APR 

178 

184 

+  6 

3.4 

MAY 

255 

260 

+  5 

2.0 

JUN 

233 

277 

-f44 

18.9 

JUL 

231 

227 

-  4 

1.7 

AUG 

241 

212 

-29 

12.0 

SEP 

236 

225 
Total 

-11 
-19 

4.7 

Average  per  month 

-  1.6 

11.3 

Below  Colusa  Basin  drain  (Table  4),  the  monthly  average  error  was 
11.3  percent,  with  a  maximum  error  of  34.3  percent,  occurring  in 
December.   The  recorded  EC  values  used  for  this  station,  which 
were  derived  from  those  published  by  the  U.  S.  Geological  Survey 
(20)  for  Knights  Landing,  were  the  monthly  averages  of  daily  sam- 
pling data  weighted  by  flow.   The  values,  therefore,  should  approx- 
imate the  true  monthly  mean  EC.   The  USGS  data  report  notes, 
however,  that  mixing  of  irrigation  return  water  with  river  water 
is  not  complete  at  the  sampling  site. 

At  Sacramento  (Table  5),  the  monthly  average  error  was  l8.4  per- 
cent, with  a  maximum  error  of  33.7  percent  for  June.   Recorded 
EC  values  for  May  through  September  were  the  monthly  averages 
derived  from  USGS  daily  sampling  data  (20)  and  should  closely 
approximate  the  flow-weighted  monthly  mean.   The  values  for 
October  through  April  are  results  of  one  day's  grab  samples  and 
could  vary  considerably  from  monthly  mean  values.   Furthermore, 
the  Sacramento  station  was  not  located  far  enough  below  the  con- 
fluence with  the  American  River  to  allow  complete  mixing.   For 
that  reason,  the  sampling  station  was  later  relocated  at  Freeport. 

Examination  of  the  actual  deviations  in  micromhos  between  the 
recorded  and  computed  values  reveals  additional  important 
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TABLE  5.   ERROR  ANALYSIS,  1955  WATER  YEAR 
Sacramento  River  at  Sacramento 


Conductivity 

in  Micromhos 

Deviation 

Month 

Recorde 

dl/ 

Computed 

Micromhos 

Percent 

OCT 

193 

147 

-47 

23.8 

NOV 

136 

174 

+38 

27.9 

DEC 

214 

199 

-15 

7.0 

JAN 

195 

152 

-43 

22.1 

FEB 

217 

162 

-55 

25.4 

MAR 

176 

151 

-25 

14.2 

APR 

131 

152 

+21 

16.0 

MAY 

153 

176 

+23 

15.0 

JUN 

175 

234 

+59 

33.7 

JUL 

209 

254 

+45 

21.5 

AUG 

224 

233 

+  9 

4.0 

SEP 

253 

228 
Total 

-25 
-15 

9.9 

Averag 

e  per  month 

-  1. 

3 

18.4 

1/  Recorded  values  for  October  through  April  are  results  of 
grab  samples  taken  on  only  one  day  each  month,  and  do  not  nec- 
essarily represent  monthly  mean  values.  Mean  monthly  flow  for 
February  was  18  percent  higher  than  flow  when  sample  was  taken. 


information.   At  Hamilton  City,  the  algebraic  total  of  monthly 
deviations  is  +66  micromhos,  an  average  of  5«5  micromhos  per 
month.   Omitting  the  nonrepresentative  sample  for  December,  the 
total  becomes  -11  micromhos,  an  average  of  only  -O.9  micromhos 
per  month. 

Below  Colusa  Basin  drain  and  at  Sacramento,  near  balances  are 
achieved  in  algebraic  totals  of  deviations  in  micromhos,  being 
-19  and  -15  micromhos,  respectively.   These  amount  to  only  -1.6 
and  -1.3  micromhos  per  month,  respectively. 

These  analyses  of  actual  deviations  show  that  the  model  computed 
very  nearly  salt-balanced  conditions  for  the  year.   Further  inves- 
tigation of  the  historic  monthly  distribution  of  salts  in  the 
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system  possibly  would  reveal  monthly  or  annual  patterns  of  salt 
distribution  which,  if  incorporated  in  the  model,  would  result 
in  more  accurate  monthly  values. 

More  importantly,  measured  EC  values  may  not  represent  a  complete 
mixture  of  all  inflowing  waters  at  some  stream  locations.  Also, 
at  those  locations  where  only  one  "grab"  sample  was  obtained  per 
month,  the  relationship  between  measured  EC  and  mean  monthly  EC 
was  uncertain. 

At  the  few  locations  where  continuous  conductivity  recorders  were 
in  operation  during  the  periods  under  consideration,  errors  in 
obtaining  mean  monthly  EC  values  were  minimized.  Reasonably  re- 
liable data  also  were  obtained  from  a  few  stations  where  samples 
had  been  collected  daily. 

The  greatest  chance  for  error  occurred  in  trying  to  account  for 
unmeasured  valley  floor  accretions.  For  those  inputs  the  sources 
were  assumed  and  values  assigned  primarily  on  the  basis  of  judg- 
ment.  For  example,  if  the  source  were  assumed  most  likely  to  be 
groundwater  inflow,  an  EC  value  corresponding  to  known  EC  of 
groundwater  in  the  area  was  assigned. 

Further  refinement  of  the  model  and  additional  evaluation  over  a 
series  of  years  of  the  more  critical  input  values,  particularly 
those  of  the  major  accretions,  undoubtedly  would  produce  more 
valid  results,  A  more  thorough  statistical  approach  toward  eval- 
uating the  relationships  between  computed  and  recorded  data  cur- 
rently is  being  Investigated.  Additional  computer  runs  with 
output  for  additional  years  would  more  materially  aid  in  the 
statistical  analysis. 

The  Department  does  plan  to  conduct  additional  verification  runs 
of  the  model  based  on  more  recent  data,  including  that  for  stations 
where  conductivity  recorders  provide  continuous  EC  records.  These 
additional  funds  could  not  be  made  within  the  time  available  for 
the  present  study  =, 

Although  additional  studies  are  desirable,  the  model  appears  to 
be  conceptually  valid.  The  I961  and  1965  runs  demonstrated  that 
a  mathematical  simulation  of  electrical  conductivity  in  the 
Sacramento  River  system  was  feasible.  The  1955  verification  run 
showed  that  the  assumptions,  coefficients  and  general  rationale 
employed  in  developing  the  model  were  sufficiently  reliable  to 
simulate  the  general  monthly  pattern  of  EC  at  selected  river  sta- 
tions and  that  approximate  values  of  mean  monthly  EC  could  be 
estimated. 

In  many  instances,  flow  and/or  EC  values  of  certain  model  inputs 
must  be  specified  primarily  on  the  basis  of  Judgment.  However, 
considerable  (and  growing)  backlog  of  data  is  available  as  guides 
to  Judgment. 
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Additional  Comments  on 
Accuracy  and  Use  of  the  "Hodel 

Output  from  the  model  consists  of  specific  values  of  flow,  EC 
and  salt  loading.  These  values  can  be  expected  to  deviate  by 
varying  amounts  from  prototype  values  where  comparisons  can  be 
made.  The  more  Important  consideration  Is  whether  or  not  the 
model  reasonably  simulates  the  prototype  system.  If  It  does, 
the  model  can  then  be  used  to  compare  alternative  systems. 

In  this  study  comparisons  were  made  In  terms  of  EC  between  two 
levels  of  cultural  development.   Results  of  the  1955  verifica- 
tion run  Indicated  that  the  model  simulated  the  prototype  system 
well  enough  to  be  used  for  estimating  the  effects  of  a  projected 
future  level  of  development  on  EC  In  the  Sacramento  River. 

The  accuracy  of  computed  results  depends  on  the  accuracy  of  In- 
put data.   In  this  study  Input  data  had  to  be  extrapolated  or 
estimated  where  measured  data  were  not  available.  Although  cer- 
tain guides  to  Judgment  were  available,  accuracy  of  those  values 
was  uncertain. 

Measured  streamflow  and  EC  values  used  in  the  model  generally 
were  accepted  as  correct,  although  they  are  subject  to  the  nor- 
mal errors  in  measurement.   These  errors,  particularly  in 
streamflow,  generally  do  not  exceed  10  percent. 
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ESTIMATES  OF  FUTURE  SALINITY  OF  WATER 
IN  THE  SACRAMENTO  RIVER  SYSTEM 


Since  the  basic  objective  of  the  study  was  to  develop  a  tool  for 
estimating  EC,  or  salinity,  under  future  levels  of  development, 
the  initial  application  of  the  model  was  placed  at  the  1990  level. 
The  selection  of  1990  was  based  on  the  availability  of  estimates 
(5)  of  anticipated  levels  of  population  growth  and  agricultural, 
municipal,  and  industrial  development  for  that  year.   The  objec- 
tive was  to  illustrate  the  extent  to  which  potential  changes  in 
any  or  all  these  areas  may  affect  the  mineral  quality  of  the 
Sacramento  River  some  20  or  more  years  in  the  future. 

Water  supply  estimates  (2)  for  1990  in  the  Valley  prepared  by 
the  Department  of  Water  Resources  and  the  U.  S.  Bureau  of  Recla- 
mation provided  a  base  hydrologic  record  and  reservoir  operation 
model  for  this  study.   Their  estimates  covered  the  33-year  period 
from  (and  including)  1922  through  195^,  which  included  the  years 
from  1928  to  1934,  the  most  critical  dry  period  ever  recorded  in 
the  Valley.   The  record  also  covered  several  normal  and  extremely 
wet  years.   The  time  allotted  to  this  study  restricted  evaluation 
of  historic  and  probable  future  EC  to  only  three  of  those  33  years: 
1931,  a  critically  dry  year;  I936,  a  relatively  normal  year;  and 
1938,  an  extremely  wet  year. 

The  1990  level  of  development,  which  had  received  the  greatest  at- 
tention in  water  supply  and  project  operation  studies  in  the  Vallej^ 
was  particularly  applicable  to  this  study  because  it  had  been  used 
in  previous  water  quality  model  runs,  thereby  making  more  input 
data  available  for  the  current  runs. 

Increased  land  use  in  the  Sacramento  Valley  between  I96I  and  1990 
is  illustrated  below. 


Land  Use 


1961 


1990 


Percent 
Increase 


Urban  acreage       l48,800       239,900        6I 
Irrigated  acreage    881,500      1,269,100        44 


Detailed  tabulations  for 
are  presented  in  Appendix 
ulation  in  the  Sacramento 
932,000  in  i960  to  2,470, 
was  used  as  the  basis  for 
waste  discharge  quantitie 
data  had  been  collected  d 
model  to  estimate  histori 
quired  before  future  cond 


these  and  other  years  used  in  this  study 
C,  "Present  and  Future  Land  Use".   Pop- 
Basin  is  expected  to  increase  (5)  from 
000  in  1990.   This  I65  percent  growth 

projections  of  municipal  and  industrial 
s  for  1990.   Because  very  little  quality 
uring  1931,  1936  and  1938,  use  of  the 
c  EC  conditions  for  those  years  was  re- 
itions  could  be  assessed.   Projected 
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1990  conditions  were  then  superimposed  upon  the  historic  conditions 
to  make  the  projected  EC  estimates. 

The  major  changes  in  the  Sacramento  Valley  water  resource  system 
between  historic  and  1990  conditions  were  the  effects  of  seven 
major  reservoirs  on  flow  into  the  Sacramento  River.   These  were 
Shasta  Lake,  Lake  Oroville,  Polsom  Lake,  Whiskeytown  Reservoir, 
Marysville  Reservoir,  New  Bullards  Bar  Reservoir,  and  Auburn  Res- 
ervoir. Consideration  also  had  to  be  given  the  system  from  the 
Trinity  River  via  Clear  Creek  and  Spring  Creek  Tunnels  and  the 
possibility  of  additional  imports  from  the  Eel  River.   For  this 
study.  Eel  River  imports  were  not  considered,  but  they  may  be 
included  in  later  computer  runs. 

Other  significant  changes  were  increased  agricultural,  municipal, 
and  industrial  development.  Agricultural  development  is  incor- 
porated in  the  accretion  analysis  of  the  study,  which  reflects  in- 
creased Irrigated  acreage  and  projected  cropping  patterns  (see 
Appendix  C).   Waste  discharges  from  municipal  and  industrial  devel- 
opment were  considered  separately  in  the  model  and  estimates  of 
future  quantities  were  based  primarily  upon  the  Department's  pop- 
ulation projections  for  use  in  water  supply  studies. 

Figures  9  through  l4  illustrate  historic  and  projected  1990  flow 
and  EC  at  Colusa  and  Freeport,  two  key  stations  on  the  Sacramento 
River.   Figure  I5  summarizes  projected  1990  results  for  the  three 
hydrologic  conditions  for  the  Sacramento  River  station  at  Freeport. 

Tables  6  through  n  are  copies  of  computer  printouts  showing  esti- 
mated monthly  flow  in  thousands  of  acre-feet,  mean  monthly  con- 
ductivity in  micromhos  and  monthly  tons  of  salt  contained  in  the 
flow  at  all  10  stations  on  the  Sacramento  River  and  on  inflow  from 
the  Feather  and  American  Rivers.   Tables  6  through  8  show  com- 
puted historic  values,  and  Tables  9  through  11  show  projected 
1990  values  for  each  of  the  three  historic  hydrologic  conditions. 
Historic  flow  values  shown  for  "Keswick  Out"  in  the  tables  were  the 
natural  flows  that  occurred  at  the  Shasta  Dam  site. 

As  indicated  earlier  under  the  heading  "Model  Verification",  the 
absolute  values  computed  by  the  model  in  some  months  may  deviate 
substantially  from  prototype  values.  However,  when  the  two  levels 
of  development  are  compared  by  means  of  computed  values  for  each 
level,  the  degree  of  difference  of  change  between  the  two  levels 
should  be  reliably  reflected  by  the  model. 

Discussion  of  Results 

This  section  compares  historic  and  projected  EC  values  for  each 
of  the  three  water  years  considered.   These  comparisons  should 
provide  some  indication  of  the  degree  of  change  which  may  be 
expected  at  the  1990  level  of  development.   The  data  are  shown  in 
Tables  6  through  11. 
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1931  Hydrology 

As  mentioned  earlier,  I93I  was  one  of  the  driest  years  on  record 
for  the  Sacramento  Valley.   The  year's  total  outflow  at  Sacramento 
was  about  4l  percent  of  that  which  occurred  in  1936,  approximately 
a  normal  water  year  (Tables  6  and  7).   Flow  values  shown  on  the 
tables  do  not  include  overflows  into  Yolo  Bypass.   During  a  dry 
year,  when  lesser  amounts  of  water  are  available  for  dilution, 
higher  than  normal  concentrations  of  dissolved  minerals  are  ex- 
pected in  river  flows. 

In  1931  maximum  computed  EC  at  Freeport  was  584  micromhos,  occur- 
ring in  August  (Table  6).   Projecting  1931  hydrologic  conditions 
to  the  year  I990  shows  a  marked  decrease  in  maximum  EC  to  362 
micromhos,  in  April  (Table  9  ).   The  August  I99O  value  dropped  to 
259  micromhos.   The  reasons  for  these  changes  are  apparent  in  view 
of  the  respective  flow  values  for  each  level  of  development.   In 
1931,  flow  in  August  at  Freeport  was  only  33,000  acre-feet;  pro- 
jected flow  for  August  1990  is  259,000  acre-feet.   On  the  other 
hand,  flow  in  April  I93I  was  492,000  acre-feet;  projected  flow  for 
April  1990  is  362,000  acre-feet. 

The  flow-weighted  yearly  average  EC  at  Freeport,  however,  was 
171  micromhos  in  I93I  and  a  projected  value  of  221  micromhos  in 
1990.   This  is  an  increase  of  about  30  percent.   However,  the 
projected  overall  mineral  quality  is  still  quite  good  and  does 
not  approach  the  historic  extremes. 

The  predicted  EC  values  for  all  stations  (Table  9  )  indicate  total 
dissolved  solids  concentrations  well  below  the  recommended  limit 
of  500  parts  per  million  TDS  for  drinking  water  set  by  the  U.  S. 
Public  Health  Service.   In  terms  of  EC,  that  limit  represents 
about  700  micromhos.   The  computed  values  are  even  farther  below 
established  criteria  for  irrigation  water.   In  a  water  supply, 
however,  the  greater  the  margin  of  dissolved  mineral  concentrations 
under  limiting  criteria  for  beneficial  uses,  the  more  suitable  is 
the  water  for  reuse.  Also,  as  mineral  concentrations  increase, 
more  water  is  needed  to  remove  salts  from  a  basin  and  maintain  a 
favorable  salt  balance, 

1936  Hydrology 

Water  year  1936  is  considered  a  normal  year  in  terms  of  Sacramento 
Valley  outflow.   In  historic  I936  maximum  EC  at  Freeport  was  esti- 
mated at  333  micromhos,  occurring  in  August  (Table  7).   The  pro- 
jected 1990  maximum  EC  at  Freeport  based  on  I936  hydrology  was 
239  micromhos,  occurring  in  June  (Table  10). 

The  flow-weighted  average  EC  for  the  year  at  Freeport  for  I936  was 
128  micromhos.   The  projected  I99O  flow-weighted  annual  average, 
based  on  I936  hydrology,  was  188  micromhos,  an  increase  of  47  per- 
cent (Tables  6  and  10).   As  indicated  in  projections  based  on  1931, 
the  overall  mineral  quality  for  199O  is  quite  good.   The  more 
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highly  regulated  river  flows  tend  to  smooth  out  quality  fluctu- 
ations and  reduce  extreme  values. 

1938  Hydrology 

In  contrast  to  1931>  water  year  I938  was  one  of  the  wettest  on 
record  in  the  Sacramento  Valley.  Not  counting  weir  spillage  into 
Yolo  Bypass,  outflow  at  Freeport  for  1938  exceeded  that  of  1936 
by  60  percent  (Tables  7  and  8).   In  historic  I938  maximum  EC  at 
Freeport  was  estimated  to  be  23I  micromhos,  occurring  in  August 
(Table  8).  The  projected  199O  maximum  EC  at  Freeport  for  1938 
hydrology  was  289  micromhos,  also  occurring  in  August  (Table  11), 

The  flow-weighted  1990  average  EC  for  the  year  at  Freeport  was 
160  micromhos,  an  increase  of  47  percent  over  the  historic  flow- 
weighted  1938  EC  of  109  micromhos.  As  expected,  the  larger  trib- 
utary inflows  reduced  I938  average  EC  of  1931  and  1936. 
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Figure  9.   HISTORIC  AND  PROJECTED  1$90  FLOW  AND  CONDUCTIVITY, 
SACRAMENTO  RIVER  AT  COLUSA,  1931  HYDROLOGY 
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Figure  10.   HISTORIC  AND  PROJECTED  1990  FLOW  AND  CONDUCTIVITY, 
SACRAMENTO  RIVER  AT  COLUSA,  1936  HYDROLOGY 
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Figure   11.      HISTORIC   AND   PROJECTED   1990   FLOW   AND  CONDUCTIVITY, 
SACRAMENTO  RIVER   AT  COLUSA,      1938  HYDROLOGY 
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Figure  12.   HISTORIC  AND  PROJECTED  1990  FLOW  AND  CONDUCTIVITY, 
SACRAMENTO  RIVER  AT  FREEPORT,  1931  HYDROLOGY 
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Figure  13.   HISTORIC  AND  PROJECTED  1990  FLOW  AND  CONDUCTIVITY, 
SACRAMENTO  RIVER  AT  FREEPORT,  1936  HYDROLOGY 
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Figure  l4 .   HISTORIC  AND  PROJECTED  1990  FLOW  AND  CONDUCTIVITY, 
SACRAMENTO  RIVER  AT  FREEPORT,  1938  HYDROLOGY 
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Figure  15.   PROJECTED  1990  FLOW  AND  CONDUCTIVITY,  SACRAMENTO 
RIVER  AT  FREEPORT,  1931,  1936,  AND  1938  HYDROLOGY 
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TABLE  6.  ESTIMATtU  1931  MONTHUY  FLO«. 


SACRAMCNTO  Rfvei)  SYSTEM 


1 

KESKlC*  OUT 

18* 

? 

»tO  HUUFF  1)0 

204 

3 

MAMILTON  CTT 

ZIU 

* 

COLUSA 

S 

Bf.L  rtUTTE  C» 

b 

»«OVt"  CBO 

r 

MLLO<  Clio 

n 

FUtMEs-NIC 

9 

VFaUN* 

10 

»MrrfIC«N-S*C 

1' 

n 

SACH4HENT0 

A14 

12 

FHteCOWT 

*22 

2*6 
2S6 

2*» 
tit 

lit 


272 

2*6 

236 

1»* 

*sr 

2TII 

216 

5*7 

23g 

li« 

S2» 

2*a 

112 

53'' 

2T3 

156 

Sk» 

1»2 

m 

5?U 

217 

111 

2'J 

122 

63 

6V* 

76 

T7» 

77» 


132 

VII2 
til 


157 
*9« 


157 

15» 

CONOUCTlvIlT  In  micmommOS 


t 

kfswick  out 

2 

HEO  rtLUFF  UO 

maxIlTON  CTT 

COLUSA 

HtL  UUTTE  CR 

AHOVE  CBO 

rtfLO*  CBO 

FFATMER-NIC 

Vt«ON« 

10  AMEHICAN-SAC 

11  SAC«»«FNTO 

12  FHf.Kt>OHT 


1*2 

1*2 
1*2 
1*2 
105 
126 


138 
137 
1«* 
1** 
lb* 
167 
167 
108 
155 


153 
155 
153 
172 

17* 


131 

1*0 
1*3 
1*5 
1** 
152 
152 
lU 
1«* 


125 
127 
132 

135 
1*0 
1*0 


130 
133 
135 
200 
211 
211 
31* 


162 
250 
250 
518 
152 
*02 


138 

1*0 

156 

152 

153 

150 

222 

300 

270 

3*6 

270 

3*6 

609 

790 

235 

1021 

5** 

687 

79 

133 

♦  02 

310 

*09 

388 

1*1 

152 
162 
10* 
265 
26S 
715 
205 
626 
100 
573 
58* 


SALT    LOAIlINft    IN    100    IONS 


1 

KFsmCK  OUT 

252 

2 

rtf-0  HL'JFF  00 

299 

3 

HAMILTON  CTT 

30' 

* 

COLUSA 

J2B 

5 

H£L  BUTTE  CH 

368 

6 

AMOVE  CBO 

365 

T 

BFlO*  CBO 

«2» 

8 

FtATHErt-NlC 

203 

9 

Vf hunA 

673 

10 

AMEhICAN-SAC 

IS 

11 

S»C"a'<ENIO 

691 

12 

FBEtHOBT 

712 

259 
310 
3*3 

351 
355 
3S5 

3d5 
202 
557 


356 

3*8 

*S6 

693 

522 

568 

835 

597 

680 

815 

616 

701 

825 

63* 

781 

9  78 

68* 

750 

991 

68* 

757 

311 

278 

*19 

1198 

1002 

no* 

*8 

56 

8* 

1261 

109O 

1188 

1282 

nil 

1202 

307 
369 
312 
*95 
576 
»05 
681 
116 
8** 


2»» 

309 
227 
181 
390 

203 


239 
292 
1T2 
22* 
338 


221 
2*2 

120 
103 


172 
193 


216 
233 

176 
172 
252 
206 
311 


TABLE  7.  ESTIMATED  1936  MONTi^LT  FLOXi 


MEAN  CONOUCTIVITY  ANO  SALT  LOaDINO  I» 
19)6  HTOnOLUSr 


sacamento  bive"  system 


FLOW  IN  innn  ACHE  FEET 
ID    station        OCT 


KFSWICK  OUT 

181 

178 

212 

805 

1062 

562 

503 

295 

202 

170 

IIKO 

HEO  HLIIFF  DO 

20« 

21* 

286 

1*26 

1667 

762 

685 

370 

233 

1H8 

177 

MAMILTnN  CTY 

22* 

2J3 

315 

17^6 

20H* 

878 

79« 

192 

1.2 

1S9 

COLUSA 

215 

2J7 

312 

l**2  . 

1713 

9ll 

788 

152 

lOl 

130 

UEL  BUTTE  CB 

23* 

2*7 

3»1 

16*6 

2116 

973 

837 

192 

1*2 

175 

AHOue  CBU 

232 

2*7 

3*1 

1**5 

1879 

89* 

762 

131 

83 

139 

RFlOW  CBO 

235 

2*7 

3*1 

1605 

2073 

90* 

762 

1*1 

103 

177 

FEATmEo-nIC 

10» 

97 

138 

599 

89* 

589 

700 

135 

131 

125 

VCn'lNA 

33j 

3<>* 

»79 

1*** 

1*36 

13*8 

1.61 

309 

27* 

337 

A-!FMIC«N.<;*C 

33 

OS 

51 

♦  29 

809 

♦  32 

627 

83 

27 

?5 

SaCJAMFNTO 

361 

391 

533 

\9U2 

1987 

17b3 

2P86 

1*»T 

930 

371 

286 

FhEEBOMT 

36* 

39* 

536 

1905 

1990 

1786 

2n«9 

1*70 

933 

37* 

289 

3')0 

CONDUCTIVITY  IN  MICBOMMOS 


1  KESWICK  OUT 

138 

139 

137 

120 

113 

122 

120 

12* 

123 

135 

138 

139 

2  REO  BLUFF  00 

158 

1*1 

137 

13* 

123 

129 

132 

125 

125 

UO 

1»7 

1*» 

3  HIVMILTON  CTY 

170 

1*5 

1*2 

1*2 

175 

130 

1** 

169 

ISO 

200 

199 

4  COLUSA 

179 

U7 

1*5 

1*6 

13* 

138 

IS? 

178 

195 

218 

206 

5  BF.L  BUTTE  CB 

200 

151 

150 

1*5 

132 

1*5 

17* 

213 

260 

301 

?25 

6  AHOVE  CHO 

200 

151 

150 

1*5 

133 

1*5 

17* 

213 

260 

301 

225, 

7  BFL'IW  CBO 

210 

151 

150 

1*6 

150 

1*9 

18* 

2*6 

298 

398 

2*6 

8  FFATmFB.nIC 

121 

125 

169 

77 

73 

77 

93 

97 

10* 

10* 

109 

9  vrpniiA 

182 

U* 

155 

127 

128 

121 

I  10 

158 

2p6 

10  ♦MF'JicAN-sie 

85 

79 

81 

57 

52 

56 

50 

55 

70 

KO 

HO 

11  SacxamknTO 

177 

139 

153 

113 

96 

106 

117 

161 

2** 

329 

213 

12  FRfEPODT 

182 

1** 

156 

11* 

97 

107 

118 

163 

2*8 

333 

?|7 

SALT  LOAOINn  IN  100  TONS 


1  KESWICK  OUT 

250 

2*7 

290 

96* 

1200 

686 

60* 

363 

273 

235 

222      i 

2  REP  HLuFF  DO 

322 

401 

191* 

2050 

983 

90* 

*6* 

326 

26^     e 

3  HAMILTON  CTY 

382 

33H 

2*.5 

26  0* 

11*3 

1057 

627 

3*6 

♦  COLUSA 

385 

3^9 

210* 

2301 

1256 

no* 

612 

296 

5  BFL  aUTTE  CH 

*67 

373 

2383 

2792 

1*15 

1182 

821 

♦  99 

39^     U 

6  ABOVE  CBO 

*63 

373 

2093 

2*96 

1301 
1351 

1076 

691 

3^0 

313     10 

7  BELOW  CBO 

*93 

373 

23*0 

3117 

lc76 

♦20 

8  FFaThER-NIC 

126 

ill 

♦  61 

65* 

♦  56 

539 

2*1 

MO 

13*      ] 

9  VFRCNA 

610 

1835 

1831 

1632 

161* 

i*28 

131* 

739     1« 

10  AMERICAN. SAC 

28 

36 

2^6 

*22 

2*0 

317 

1«1 

58 

20 

n  SACBAMFNTO 

6*0 

5^5 

2158 

1,0, 

1*93 

I95I 

1718 

1502 

,06 

12  FBFEPOOT 

661 

566 

2179 

l»3o 

191* 

1972 

1739 

1523 

927 

9*2 

781     IS 

48 


TABLC  8.  CtrlMATtO  !«]■  MONTHCr  fkOVi 


FLOW  IN  Idle  *C«e  f«T 
10    ST»TIOH        »CI 


Ke^mCK  ojT 

H5 

T»» 

Hto  RLi/rr  uo 

f3* 

n»T 

MiMILTON  CTT 

ZSl 

ij^i 

COLOS* 

itr 

11** 

BfL  UDTTE  C« 

Mil 

w* 

»Mnvr  CMO 

»5A 

loat 

HULO*  C«0 

2*« 

1011 

FeATHCD-NtC 

li* 

»i<i 

VOUNt 

*ll 

1012 

It 

IflEMlCiN-Ste 

]« 

«• 

II 

S«C«»MFNTO 

»»<! 

Utl 

12 

FBeCPOBT 

♦SJ 

l*k* 

zieo 

■«C*N  COMOUCTIVITV  talO  tALT  LOtOIN*  IN  l*C**MCNTO  KIVCK  STtTCM 


!»♦» 

ir*s 

21*0 
ITOO 
1«IZ 

io*z 

»«« 

IaS3 

,"« 
ZJTS 

Z)T( 


S«T 

1»3| 

!••* 

»10 

t9«« 

ilfo 

)!»» 

l»«» 

20  04 

l»" 

?5*» 

imo 

I»»l 

ij»« 

2*21 

»i3 

1S2» 

1946 

11131 

P* 

S7« 

list 

2nin 

Jl»l 

2021 

*ll* 

*T* 

101 

2*T 

223 

*•*• 

w*> 

9*S 

1** 

2»» 

2»T 

l*l»» 

!»»♦ 

T23 

J3» 

22* 

23* 

ITMt 

H»» 

*2* 

312 

201 

20f 

1*1** 

!?♦• 

y»* 

3*0 

t»* 

2*» 

l»»*» 

1300 

t20 

272 

iro 

213 

i2Jjr 

1300 

*2l 

302 

i»i 

231 

133*» 

22** 

»0«2 

2I» 

m 

12*3* 

130* 

Uao 

SI* 

24  ♦ 

2fT 

U2»0 

loo« 

its 

12* 

34 

**2* 

3i«T 

2<>it 

*3l 

2TT 

2*3 

l«9l» 

S8»0 

2»»l 

*3* 

tut 

ta* 

lt»»» 

CONDUCTIVITY  IN  MICROlJHO* 

1  KtSWICK  OJT  I3T 

2  BFn  HtllFF  00  1*4 

3  HJMlLTnN  CTY  2SH 
•  *  COLJS*  2*^ 

5  HEL  HUTTE  CB  300 

6  »40ve  CUD  300 

7  HFLOW  CBO  1*7 
«  Ft»TMFo.N|C  1»* 
9  VFHON*  2»» 

10  *MeHlC>N>S*C  B* 

11  S^CSAMFNTO  272 

12  FuEKpnuT  271 


10» 
1*0 
144 

14lt 
ltl7 
lt7 


114 
13S 
130 
132 
13* 
13* 
1*» 


ll» 
13S 
133 

lit 
1J7 
137 

1*2 
103 


10» 
123 
121 
130 
127 
124 


II* 
12* 
123 

130 

130 
131 
134 


11* 
120 
121 
12» 
I3S 
13» 
13* 


114 
107 

lis 

117 
117 
117 


114 

114 

no 

114 
127 
127 
128 


127 
133 

12* 
13S 

140 
140 
210 


130 
1J7 
1*4 
177 
21S 
21S 
240 
130 
234 


132 

134 
1*1 
ISl 
IT* 
17* 
220 

uo 

143 
43 

1H5 

HI 


S*LT  LOJniNR  IN  100  TONS 


1  KFSJiCx  OUT 

2  PEl)  NLhFF  do 

3  HAMIUTDN  CTY 
*  COL'JS* 

5  B»L  HUTTE  C« 

6  »nnvE  CRO 

7  MFLOi.  cno 

g  FE«TMSu-NlC 

9  VFHdN* 

10  »ME"ICAN-S»C 

11  "HCWAWNTO 

12  FWEFMORT 


251 

»*« 
42? 
773 
T6T 


7BS 
143* 
22** 
1950 
22J'» 
IHMS 
22»? 

SIS 

19x0 

hi 


12*4 
2^17 
321» 
28*0 
3119 
278» 
3135 
11*3 
2279 
25* 
2221 
22*2 


49V 
1254 
I49S 
1513 
17*« 
1412 
1499 

5*0 
2417 

100 
27*} 
27<»* 


1404 
31»2 
34»0 
2710 
32«* 
237* 
3092 
ll"! 
2021 
31* 
1912 
1933 


2154 

3897 
■*720 
3634 
*0>15 
302IS 
3*37 
1*97 
1154 
*20 


Uil 
2101 
2S*2 

2196 
2**7 
1*3* 
1*3* 
l«1S 
1*** 
36S 
1H31 
1*52 


1063 

i**o 

1702 
1427 
?0S2 
1526 
I  526 
137* 
7S99 


552 

707 

796 

73* 

90* 

790 

795 

T»« 

196* 

298 

21*7 

2168 


342 

*60 
♦  21 
«21 
590 
*36 
633 
201 
831 


321 
369 
382 
36* 
551 
166 
Sb6 
105 
604 


29* 

1087* 

3*3 

14291 

32» 

22S79 

310 

19123 

*42 

2228 1 

370 

1737S 

523 

2019r 

12* 

90*9 

*9S 

l9*S9 

2*60 

57* 

2ISU 

sts 

2176* 

TABLE  9.  tSTIMATtO  1999  MONTNLY  fLOKi 


1 

KESWICK  OUT 

2*3 

BED  BLUFF  00 

259 

H«NILTnN  CTY 

2*0 

COLUS* 

BFL  bUTTe  CR 

«HOve  CBO 

BELOW  CBO 

FE»TMEH-NIC 

VFBON* 

in 

»MFBie«N.S»C 

u 

S«CB»MENTO 

12 

FHEFPOBT 

5*S 

79 

97) 

S44 

49  g 

160 
293 
i»» 
S'l 

347 
3*4 
379 
154 

56* 


171 

404 

37* 

243 

663 

**4 

30* 

m 

277 

2«« 

2»* 

i«i 

736 

341 

256 

6»7 

2*0 

257 

644 

38* 

124 

151 

45 

4*1 

817 

385 

*7 

»* 

30 

552 

8*7 

3l6 

547 

894 

331 

285 
3SS 

237 
2^5 

3*9 
2*8 

30* 
122 

*50 


243 

3*9 
2*0 
244 
3*2 
234 
317 


*77 
*92 


33* 
384 
234 
2*» 
3*« 
255 
32* 
114 
57* 
153 
59* 
41* 


334 
347 
26* 
141 
29* 
2*8 
240 
251 
560 
139 
657 
672 


142* 

*322 

3451 
)44'> 

*2*6 
3556 

4010 
191)i 
4073 

1203 
7066 

72*0 

CONDUCTIVITY  IN  MICBOMMOS 


1  KfSmCK  OUT 

12* 

130 

131 

127 

115 

2  BED  BLUFF  00 

135 

160 

13* 

143 

3  HAMILTON  CTY 

162 

147 

1*7 

272 

*  COLUS* 

17* 

193 

151 

25* 

5  BEL  BUTTE  CB 

189 

1S9 

1*9 

302 

6  «HOVE  CBO 

189 

192 

1*9 

302 

7  BELOW  CBO 

219 

250 

159 

*51 

e  FEATHCB-NIC 

114 

120 

99 

79 

9  VEBON* 

17S 

214 

1*4 

10  «UFaic«N-$«C 

11  S4CB4NENT0 

*3 

44 

96 

88 

141 

134 

212 

150 

284 

12  FBEEPOBT 

174 

1** 

223 

15? 

2*9 

t4LT  LOAOlNS  In  10*  TONS 


1  KESWICK  OUT 

114 

1*4 

20* 

224 

•«« 

«** 

14] 

IM 

4*4 

2  BCD  BLUFF  00 

150 

4*4 

1*4 

»12 

444 

»1* 

147 

545 

1  HiNRTON  CTY 

4*1 

47* 

1074 

«47 

481 

4(2 

702 

*  eOLUS> 

4** 

*7* 

10*3 

491 

715 

477 

712 

S  BEL  BUTTt  CH 

4*1 

**-« 

1»»4 

902 

i447 

10)2 

1114 

424 

*U 

97* 

400 

754 

718 
14l4 

41* 

7  BFLOW  CBO 

4*7 

*20 

1059 

1*40 

1107 

1*8* 

4  FE4THEB-NIC 

179 

11* 

150 
12t* 

SI 

-1»T 

*• 

221 

9  VEBON* 

*43 

121* 

711 

1192 

ISII 

1454 

10  4HrB|C*N-$tC 

44 

50 

<2 

12 

4* 

1)4 

1*3 

11  StCBlMENTO 

•  71 

1147 

«*4 

1111 

10*1 

III* 

1344 

172* 

12  FBEFP09T 

474 

1*24 

1472 

fl41 

1*27 

11** 

1421 

144* 

182* 

1151 

li*73 

49 


TABLE  10.  CJTIHATIO  l»»»  MONTHLY  FLOWt 


MEAN  CONDUCTIVITY  kNO  SALT  UOAOINO  IN  SACKAMCNTO  KIVER  SYSTCH 
1*1*  MYOItOLOtY 


FLO*  IN  1000  ACnt  flET 
10    STATION        OCT 


1  KESWICK  OUT 

33Z 

363 

160 

2  RED  Buurr  ou 

36* 

ASS 

51* 

3  HAMILTON  CTT 

2»0 

51» 

TOO 

♦  COLUSA 

300 

501 

872 

5  8EL  auTTE  CM 

30T 

SIO 

«»3 

6  ABOVE  C90 

2(9 

412 

97a 

7  BELOW  C80 

3*S 

♦  15 

1077 

S  FEATHER. NIC 

222 

IBS 

♦  67 

9  VERONA 

567 

601 

lie* 

1392 

10  AMERICAN-SAC 

123 

1^0 

11  SACRAMENTO 

630 

7*3 

1077 

1728 

IZ  FREEPORT 

6AS 

7S« 

1092 

1743 

49^ 
633 
393 
3?7 
379 
250 
♦  15 
3^5 
781 
\*D 
755 
770 


5^1 

657 

440 

357 

417 

271 

452 

245 

652 

259 

68t 

1048 

699 

1063 

1005 

358 
253 
258 
282 

242 

280 
302 
515 
140 
588 
603 


CONDUCTIVITv  IN  HICRQHHOS 


KESWICK  OUT 

HEO  BLUFF  DO 

HAMILTON  CTT 

COLUSA 

BEL  BUTTE  C* 

ABOVE  CSD 

BELOW  CBD 

FEATMER-NIC 

VERONA 

10 

AMERICAN>$Ab 

11 

SACRAMENTO 

12 

FREEPORT 

SALT  LOAOINO  IN  100  TONS 


US 
121 

127 
134 
133 

133 
175 
115 
156 


121 
133 
138 
143 
143 
143 
146 
118 
137 


119 
163 
163 
171 
178 
178 
198 
115 
183 


114 
129 
133 

143 
150 
150 
171 
102 
148 
83 
145 
148 


118 

116 

13^ 

132 

Ul 

140 

162 

159 

157 

160 

157 

160 

186 

299 

93 

69 

U9 

203 

73 

77 

164 

215 

168 

223 

123 
127 
156 
171 
223 
223 
333 


118 

123 

163 

177 

243 

243 

346 

122 

301 

63 

229 

224 

239 

231 

117 
122 
14* 
U3 
198 
201 
275 
129 
268 


llT 
126 
163 
175 
194 
19* 
247 
110 
172 


1  KESWICK  OUT 

384 

395 

191 

929 

336 

282 

609 

639 

1049 

1013 

66** 

2  RED  BLUFF  OU 

463 

493 

837 

1521 

582 

546 

802 

905 

1144 

1072 

9351 

1270 

1979 

653 

512 

612 

716 

854 

804 

9368 

4  COLUSA 

438 

467 

1490 

2277 

628 

508 

559 

632 

871 

787 

9827 

446 

1678 

2646 

662 

563 

846 

1013 

1308 

1123 

120»S 

1742 

2646 

599 

381 

558 

658 

1093 

956 

10499 

614 

2134 

3377 

746 

1206 

1363 

1563 

2003 

174* 

16724 

247 

198 

539 

1584 

736 

222 

284 

299 

240 

273 

5177 

902 

810 

2174 

3318 

2079 

1*76 

1884 

1961 

2388 

2067 

20789 

64 

55 

116 

295 

117 

62 

99 

164 

249 

249 

1671 

11  SACRAMENTO 

1016 

908 

1039 

1942 

3300 

2927 

1940 

1702 

1565 

2351 

2096 

2167* 

1121 

1013 

1144 

2047 

3405 

2932 

2045 

1807 

1670 

2456 

2201 

1094 

2293* 

TABLE  11.  ESTIMATED  1990  MONTHLY  FLOW, 


MEAN  CONDUCTIVITY  AND  SALT  L0ADIN9  IN  SACRAMENTO  RIVER  SYSTEM 
1938  MYDROLOSY 


FLOW  In  1000  ACRE  FEET 
10    STATION        OCT 


1  KESWICK  OUT 

297 

232 

221 

160 

2010 

987 

1008 

957 

2  RED  HLUFF  00 

355 

636 

89* 

434 

2977 

2021 

1442 

1314 

3  HAMILTON  CTY 

300 

767 

1323 

639 

3357 

2573 

1627 

1427 

4  COLUSA 

311 

709 

1294 

599 

36,3 

2953 

1421 

1302 

5  BEL  BUTTE  CR 

318 

822 

1419 

677 

3810 

3026 

1491 

1403 

6  ABOVE  CBD 

311 

7^9 

1390 

636 

378? 

2995 

1395 

1295 

7  BELOW  C80 

395 

904 

1515 

723 

.   4011 

3119 

1553 

1437 

8  FEATHER-NIC 

225 

22^ 

1190 

530 

1239 

2099 

1384 

1467 

9  VFROnA 

AZ" 

759 

1747 

1303 

37»5 

2777 

1433 

2232 

10  AMERICAN-SAC 

86 

49 

129 

151 

300 

699 

261 

418 

u  sacramf;nto 

711 

924 

1393 

1344 

3791 

30e2 

1933 

2402 

12  FBEFPOOT 

726 

839 

1408 

1360 

3806 

3097 

1948 

2417 

496 

641 

659 

302 

619 

693 

669 

3?9 

415 

427 

385 

272 

268 

369 

334 

282 

352 

437 

388 

299 

138 

249 

249 

244 

314 

424 

403 

2A4 

911 

369 

363 

530 

1272 

636 

799 

813 

522 

301 

184 

146 

1751 

9)16 

937 

971 

1766 

1001 

952 

8*6 

CONDUCTIVITY  IN  MICROMHOS 


I  kfswick  out 

116 

117 

122 

121 

li4 

118 

116 

123 

UT 

121 

IIT 

llT 

2  RFO  BLUFF  DO 

125 

12* 

125 

169 

in. 

107 

128 

129 

126 

130 

125 

129 

3  HAMILTON  CTY 

1*7 

172 

132 

195 

116 

116 

130 

136 

189 

216 

206 

197 

4  COLUSA 

160 

176 

137 

193 

121 

124 

140 

145 

204 

223 

214 

207 

5  BEL  BUTTE  CR 

159 

179 

141 

199 

1?5 

126 

136 

159 

311 

322 

302 

2?6 

6  ABOVE  CBO 

159 

179 

141 

199 

l?^ 

126 

136 

159 

311 

322 

302 

276 

7  BELOW  CBO 

192 

214 

159 

216 

138 

137 

170 

195 

534 

508 

465 

776 

B  FEATHFn-NIC 

131 

166 

116 

118 

104 

90 

93 

88 

92 

112 

117 

109 

9  VFRONA 

170 

204 

144 

182 

132 

122 

134 

150 

233 

357 

328 

1*7 

10  AMERICAN-SAC 

69 

94 

78 

93 

70 

6n 

67 

61 

65 

73 

87 

11  SacRiMFNTo 

161 

202 

140 

173 

130 

112 

127 

137 

186 

278 

282 

isa 

12  FREEPORT 

172 

211 

146 

179 

133 

115 

132 

141 

190 

284 

289 

1*7 

SALT  LOAOINO  IN  100  TONS 


1  KFSWICK  OUT 

272 

269 

193 

2289 

1046 

1165 

1177 

581 

776 

769 

352 

2  RFO  HLirFF  00 

789 

1120 

733 

3238 

2154 

1850 

1699 

778 

900 

835 

473 

3  HAMILTON  CTY 

131* 

1749 

1181 

3902 

2*84 

2115 

1»36 

786 

922 

793 

537 

4  COLUSA 

1391 

1774 

1152 

4387 

3544 

1992 

1*97 

547 

825 

714 

5"5 

5  BFL  BUTTE  CR 

1474 

199* 

1350 

4766 

3823 

2034 

2228 

1094 

1406 

1173 

676 

6  ABOVE  CBO 

1411 

19S8 

1265 

4732 

3TT5 

1903 

2057 

429 

901 

753 

7S4 

7  BELOW  CBO 

1932 

24o4 

,l62 

55,8 

42*9 

2*33 

7900 

1675 

2, si 

1973 

9  FEATMER-NIC 

373 

1375 

627 

1299 

1899 

1293 

1299 

8*3 

415 

578 

9  VFHONA 

1053 

1547 

2S19 

2375 

5007 

3378 

2591 

3339 

2966 

2982 

2622 

1359 

10  AMERICAN-SAC 

46 

87 

117 

249 

479 

179 

281 

320 

196 

135 

127 

11  SACRAMENTO 

1144 

1668 

1950 

2329 

49I9 

3447 

2463 

ijoo 

3758 

2737 

2361 

1373 

12  FREEPORT 

1251 

1773 

2055 

2434 

5044 

3552 

2569 

3405 

3363 

2*42 

2466 

1478 

32233 
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SENSITIVITY  ANALYSES 


A  final  phase  of  this  study  was  an  additional  series  of  computer 
runs  performed  to  test  the  model's  degree  of  response,  or  sensi- 
tivity, to  varying  magnitudes  of  certain  input  data.   These  addi- 
tional runs  were  expected  to  provide   1)  a  guide  to  those  areas 
where  further  refinement  of  the  model  might  be  Justified  and 
2)  an  indication  of  the  probable  range  of  computed  EC  values  in 
the  river.   The  first  step  was  to  determine  how  sensitive  the 
model  was  to  those  variables  which  have  the  widest  range  of  prob- 
able values.   As  brought  out  by  this  study,  those  variables  are 
the  quantities  and  qualities  of  valley  floor  accretions.  A 
second  step  was  to  test  the  sensitivity  of  the  model  to  munici- 
pal and  industrial  waste  discharges. 

Valley  Floor  Accretions      ,   ., 

Historically,  the  probable  quantities  of  valley  floor  accretions, 
although  widely  variable,  appeared  more  reliable  in  their  deter- 
minations than  did  the  probable  qualities  (EC's).  Accretion  quan- 
tities could  be  related  to  probable  return  flows  from  precipitation 
and  applied  irrigation  water  in  each  hydrographic   area.  Although 
accretion  EC  could  be  related  to  the  estimated  salt  input  of  the 
source  waters,  there  was  no  firm  basis  for  estimating  how  much  of 
these  salts  would  be  returned  to  the  stream  during  any  given  month. 

In  the  first  sensitivity  runs,  therefore,  accretion  EC  values  were 
varied  with  respect  to  EC  values  used  in  the  original  computer  runs. 
These  included  runs  for  the  199O  level  of  development  under  each  of 
the  three  hydrologic  conditions--1931,  1936  and  1938.   The  previous 
EC's  of  all  major  accretions  in  the  system,  those  occurring  in 
Reaches  B,  D,  F,  G  and  H  (Figure  1),  were  multiplied  by  factors  (C) 
of  0.5,  0.9,  1.1  and  1.5  for  the  three  hydrologic  conditions. 

Figures  16  through  21  show  the  monthly  results  at  Colusa  and  Free- 
port  for  each  year  with  C  =  O.9,  1.0  and  1.1.   Figure  17,  in  addi- 
tion, shows  values  for  C  =  O.5  and  1.5.   The  0.5  and  I.5  C  factors 
were  used,  not  to  produce  a  probable  range  of  values,  but  to  show 
the  overall  response  of  EC  in  the  river  to  accretion  EC's. 

The  results  of  the  sensitivity  runs  with  C  =  0.9  and  1.1  could  be 
viewed  as  a  possible  range  of  predicted  EC  at  the  river  stations. 
This  range  would  be  based  on  the  assumption  that  each  accretion  EC 
value  as  estimated  for  the  original  runs  had  a  probable  error  of 
+10  percent.   A  more  statistically  reliable  range  could  be  obtained 
by  reexamining  each  Individual  accretion  EC  value  as  to  its  probable 
range  of  accuracy  and  then  running  the  model,  using  first  the  low 
range  and  then  the  high  range  for  each  EC.   Undoubtedly  some  of  the 
original  values  are  considerably  more  accurate  than  others  and  a 
blanket  10  percent  probable  error  for  each  EC  value  is  not  realistic, 
However,  a  lack  of  time  during  this  study  prevented  a  more  detailed 
evaluation. 
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The  figures  clearly  display  the  significance  of  valley  floor  ac- 
cretions during  the  irrigation  season,  which  is  evidenced  by  the 
wider  spread  between  each  month's  values  for  April  through  Septem- 
ber.  The  sharp  response  of  the  river  EC,  with  accretion  EC  values 
varying  +50  percent  (Figure  17)*  shows  that  the  quality  of  the 
river  can  be  profoundly  influenced  by  the  quality  of  the  accretions. 

The  results  of  these  runs  justify  additional  study  to  more  accurately 
evaluate  both  quantity  and  quality  of  valley  floor  accretions. 

Municipal  and  Industrial  Waste  Discharges 

Sensitivity  runs  also  were  made  to  test  the  effects  on  river 
quality  of  municipal  and  industrial  waste  discharges.  In  these 
runs,  all  1990  flow  quantities  of  M&I  wastes  for  1931  and  I936 
hydrology  as  originally  used  were  doubled.  EC  values  were  held 
at  the  constant  700  micromhos,  a  value  based  upon  recorded  data 
for  the  City  of  Sacramento  sewage  treatment  plant  that  probably 
will  not  change  significantly. 

With  the  original  M&I  waste  flows  doubled,  the  annual  flow- 
weighted  average  EC  for  1931  hydrology  in  the  Sacramento  River  at 
Freeport  increased  from  221  to  237  micromhos,  or  7.2  percent.   For 
1936  hydrology,  the  increase  was  from  188  to  193  micromhos,  or  2.7 
percent.   Results  for  individual  months  are  shown  in  Table  12. 


TABLE  12.  EFFECTS  OF  DOUBLING  MUNICIPAL 

AND  INDUSTRIAL  WASTE  FLOWS  ON  1990  EC, 

SACRAMENTO  RIVER  AT  FREEPORT 


(EC  in  micromhos  at  25 

°  C) 

Using  Original 
M&I  Waste  Flows 

Using  2  X  Original 
M&I  Waste  Flows 

Month 

1931       1936 

1931       1936 

OCT 
NOV 
DEC 
JAN 
FEB 
MAR 
APR 
MAY 
JUN 
JUL 
AUG 
SEP 


176 

174 

167 

175 

149 

151 

223 

187 

195 

148 

159 

168 

362 

223 

317 

235 

299 

239 

298 

231 

259 

219 

171 

181 

196 

182 

186 

183 

165 

158 

238 

192 

213 

150 

172 

171 

381 

228 

332 

240 

315 

244 

311 

235 

273 

223 

187 

189 
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Figure   l6.      SENSITIVITY   ANALYSIS,    SACRAMENTO  RIVER    AT  COLUSA, 
1990  LEVEL,    1931   HYDROLOGY 
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OCT       NOV       DEC        JAN       FEB       MAR       APR       MAY        JUN        JUL        AUG        SEP 
500 


Figure   17.      SENSITIVITY   ANALYSIS,    SACRAMENTO  RIVER    AT  PREEPORT, 
1990  LEVEL,    1931    HYDROLOGY 
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Figure  18.   SENSITIVITY  ANALYSIS,  SACRAMENTO  RIVER  AT  COLUSA, 
1990  LEVEL,  1936  HYDROLOGY  ,  ,  ■  ,   ,.   .  ^ 
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Figure  19.      SENSITIVITY  ANALYSIS,    SACRAMENTO  RIVER  AT  FREEPORT, 
1990  LEVEL,    1936   HYDROLOGY 
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Figure   20.      SENSITIVITY   ANALYSIS,    SACRAMENTO  RIVER   AT   COLUSA, 
1990  LEVEL,    1938  HYDROLOGY 
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Figure    21.      SENSITIVITY   ANALYSIS,    SACRAMENTO  RIVER   AT  FREEPORT, 
1990  LEVEL,    1938   HYDROLOGY 
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APPENDIXES 


APPENDIX  A 

ANNOTATED  BIBLIOGRAPHY 

Ref . 
No. 

(1)  California  State  Department  of  Water  Resources.   "Accretion 

Study  Phase  of  the  July  I966  DWR  Water  Rights  Operation 
Study."   Office  Report.   April  I967. 

(2)   .      "Central  Valley  Project-State  Water  Project  System 

Operations  Study  No.  90-1. 01 -69."  DWR  Files.   December 
1969. 

This  publication  provided  1990  Inflow  quantities. 

(3)   .      "Explanation  of  Input  Data,  Joint  DWR-USBR  Central 

Valley  Depletion  Study  of  July  I966."  Office  Report. 
September  1965. 

(4)  .   "Ground  Water  Conditions  in  Central  and  Northern 
CJalifornia."  Bulletin  No.  77  series. 

(5)   .   "Implementation  of  the  California  Water  Plan." 

bulletin  No.  I6O-66.   March  I966. 

(6)   .   "Reclamation  of  Water  from  Sewage  and  Industrial 

Pastes  in  California,  July  1,  1955- June  30,  1962." 
Bulletin  No.  68-62.   October  I963. 

Basic  source  of  municipal  and  industrial  waste  discharge 
data. 

(7)   .   "Sacramento  Valley  Seepage  Investigation."  Bulletin 

Uo.  125.   August  1967. 

(8)   .      "Supplemental  Information  for  DWR  Exhibit  No.  80, 

Coordinated  State  Water  Project-Central  Valley  Project 
Operation  Study."   August  1966. 

.   "Upper  Eel  River  Development  Interim  Report,  Water 

Quality."  Office  Report.  August  I966. 

.   "Water  Rights  and  Estimated  Entitlements  to  the  Flow 

of  the  Feather  River."  Bulletin  No.  l40.   August  1965. 

United  States  Department  of  the  Interior,  Bureau  of 
Reclamation;  California  State  Department  of  Water 
Resources;  Sacramento  River  and  Delta  Water  Association. 
"Hydrology  Supplement  to  Report  on  1956  Cooperative 
Study  Program,  Water  Use  and  Water  Rights  Along 
Sacramento  River  and  in  Sacramento-San  Joaquin  Delta." 
March  1958. 
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(12)   .  "1957  Joint  Hydrology  Study,  Sacramento  River  and 

?an  Joaquin  Delta."  July  195o. 

(13)  United  States  Department  of  the  Interior,  Geological  Survey. 

'Chemical  Composition  of  Snow  in  Northern  Sierra  Nevada 
and  Other  Areas."  Water-Supply  Paper  1535- J.  1964. 

(14)  University  of  California,  Agricultural  Extension  Service. 

"West  Side  Sacramento  Valley  Projected  Subsurface 
Drainage  Investigation,  Plows  and  Quality."  Unpublished 
Report.   April  1970. 

(15)  Water  Resources  Engineers,  Inc.;  California  State  Department 

of  Water  Resources.   "Develop  Preliminary  Structure  of 
Hydrologic-Water  Quality  Model."  San  Francisco  Bay-Delta 
Water  Quality  Control  Program.   Task  Order  No.  III-l. 
March  19 67. 

The  following  publications  were  basic  sources  of  streamflow 
and  water  quality  data. 

(16)  California  State  Department  of  Water  Resources.   "Hydrologic 

Data",  Appendix  B:   "Surface  Water  Plow",  Appendix  D: 
"Surface  Water  Quality",  Appendix  E:   "Ground  Water 
Quality".  Bulletin  No.  I30  Series. 

(17)   .      "Quality  of  Ground  Waters  in  California."  Bulletin 

Ro.  66  Series. 

(18)   .      "Quality  of  Surface  Waters  in  California."  Bulletin 

Flo.  65  Series. 

(19)   .   "Surface  Water  Plow."  Bulletin  No.  23  Series. 

(20)  United  States  Department  of  the  Interior,  Geological  Survey. 

"Quality  of  Surface  Waters  of  the  United  States  1955." 
Parts  9-14.   Water  Supply  Paper  l403. 

(21)   .   "Surface  Water  Records  of  California."  Series. 

(22)   .      "Surface  Water  Supply  of  the  United  States."  Water- 

S"upply  Paper  1715>  Part  II.   Series. 

The  following  publications  were  of  especial  value  in  pre- 
paring this  report. 

(23)  California  State  Department  of  Water  Resources.   "A  Machine 

Program  to  Estimate  Monthly  Consumptive  Use  and  Water 
Requirements."  Office  Report.   November  I967. 

This  program  is  commonly  referred  to  as  the  C.U.-2 
Program. 
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(24)  .   "Estimated  1990  Mineral  Quality  of  Water  In  the 

^      gacramento  River  Based  on  Dry  Period  (1931-3^)  Hydrology." 

Office  Report.  November  196? . 

(25)   .   "Predicted  1990  Mineral  Quality  of  Water  In  the 

^      Sacramento  River  with  Einphasls  on  the  Effects  of  Agri- 
cultural Drainage."  Office  Report.   November  1966. 

(26)  .   "Sacramento  River  Water  Pollution  Survey."  Bulletin 
^      JJo.  ill.   Appendix  A,  "Hydrography,  Hydrology  and  Water 

Utilization  and  Appendix  B,  "Water  Quality".   August 
1962. 

(27)  .   "Water  Quality  Management  of  the  State  Water  Project, 
^      Central  Valley  Portion,  Part  I:   Sacramento  Valley. 

Office  Report.  June  1965. 
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APPENDIX  B 
RESERVOIR  COMPUTATIONAL  METHODS 


Two  basic  approaches  were  employed  by  the  I965  model  to  predict 
the  mineral  content  of  reservoir  outflow.   Earlier  in  this  report 
computational  methods  were  expressed  verbally  but  mathematical 
equations  were  omitted  for  the  sake  of  brevity.   This  appendix 
presents  the  mathematical  equations  for  each  approach.   These 
are  identified  as  Method  1  and  Method  2. 

■   '   '    Method  1  -  :. 

The  quality  of  reservoir  releases  was  assumed  equal  to  the  qual- 
ity of  water  in  storage,  determined  by  successively  mixing  the 
flow-weighted  average  quality  of  inflow  for  each  month  with  the 
average  quality  of  water  in  storage  for  that  month,  and  then  al- 
lowing for  a  one-month  flow-through  time  in  the  reservoir.  The 
inflow  quality  is  computed  by  the  C/Q  equation  (see  "The  I965 
Model"),  and  the  initial  quality  of  water  in  storage  is  assumed 
equal  to  the  quality  of  inflow  for  the  month  preceding  the  initial 
month  of  consideration. 

In  Method  1,  a  weighted  average  is  obtained  of  the  quality  of  the 
inflow,  storage,  and  outflow.  The  quantity  of  the  total  inflow 
into  the  reservoir  may  be  calculated  from  the  following  equation. 

I  =  A  S  +  E  +  D  +  0  .     Eq .  (B 1 ) 

where: 

I  =  average  monthly  inflow  in  cfs; 
AS  =  change  in  end-of -month  storage  in  cfs  per  month; 
E  =  average  monthly  evaporation  rate  in  cfs; 

D  ■=  average  monthly  diversions  out  of  the  reservoir 

in  cfs;  and 
0  =  average  monthly  outflow  in  cfs. 

The  quality  of  the  total  inflow  may  be  assumed  equal  to  the 
weighted  average  of  the  measured  (or  calculated)  quality  of 
all  the  tributaries. 
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J  =  1  Eq.  (B2) 


Cm 


n 


I 


^J 


J  =  1 


where: 


m 


=  specific  conductance  (EC)  of  the  total  inflow 
for  any  month  under  consideration; 


q  =  average  monthly  discharge  of  a  tributary  in  cfs; 
c  =  monthly  EC  from  C/Q  equation; 
n  =  number  of  tributaries;  and 
j  =  any  tributary 

The  data  obtained  from  Equation  (B2)  are  used  as  input  data  for 
Equation  (B3)j  where  it  is  assumed  that  the  quality  of  outflow 
for  a  given  month  is  equal  to  the  quality  of  the  water  in  storage 
for  the  previous  month;  that  is,  the  outflow  for  one  month  occur- 
red before  the  inflow  for  that  month  affected  the  quality  of  the 
water  in  storage. 

The  equation  is : 

C(l)jjj  =  C(l)m-1  (Sm-l-Om-Dm)+CmIm        Eq.  (B3) 

^m 

where: 

C(l)  =  monthly  specific  conductance  obtained  from 
Method  1; 

S  =  end-of -month  storage  in  terms  of  flow  in  cfs 
required  to  produce  the  given  storage  in  one 
month ; 

0  =  mean  monthly  outflow  in  cfs; 

D  =  mean  monthly  diversions  in  cfs; 

1  =  mean  monthly  inflow  from  Equation  (b1); 

C  =  average  specific  conductance  of  the  total 
Inflow  from  Equation  (b2); 

m  =  month  of  consideration;  and 

m-1  =  month  preceding  m. 
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In  Equation  (b3)  the  original  C(l)^_]^  value  may  be  any  reasonable, 
easily  obtained  value,  since  any  error  caused  by  an  incorrect 
initial  value  diminishes  progressively  in  calculations  for  following 
months. 

Method  2 

The  quality  of  reservoir  releases  was  assumed  equal  to  the  quality 
of  water  in  storage,  determined  by  the  running  average  of  mean 
monthly  inflow  quality,  weighted  according  to  flow,  over  a  period 
of  time  in  which  the  total  volume  of  outflow  during  the  period 
equaled  or  exceeded  the  volume  of  storage  at  the  end  of  the  period. 

In  Method  2  a  weighted  average  of  the  inflow  is  averaged  over  a 
time  period  determined  by  comparing  the  reservoir  content  with  the 
reservoir  outflow.   The  period  may  be  determined  for  each  month  by 
adding  the  total  outflow  (diversions,  evaporation,  and  reservoir 
releases)  for  the  month  under  consideration  to  the  total  outflow 
for  the  proper  number  of  consecutive  months  which  immediately  pre- 
cedes  the  month  under  consideration.   The  number  of  months  com- 
prising the  time  period  should  have  been  determined  when  the  total 
volume  of  the  monthly  outflows  equals,  or  just  exceeds,  the  volume 
of  water  in  storage  at  the  end  of  the  month  under  consideration. 
This  procedure  is  expressed  in  the  following  equation. 

:  m  ■:i:x.'',r:.       ■'.■'■ 

Sm  -  y^  (0  +  E  +  D)^   .„.,  .   .......    Eq.  (b^^) 


i  =  X 


where: 


S  =  erid-of-month  storage  in  terms  of  the  flow  in  cfs 

required  to  produce  the  given  storage  in  one  month; 

0  =  mean  monthly  outflow  in  cfs; 

E  =  mean  monthly  evaporation  rate  in  cfs; 

D  =  mean  monthly  diversions  in  cfs; 

i  =  any  month  during  the  period  of  consideration; 

m  =  calendar  number  of  the  month  under  consideration;  and 

X  =  calendar  number  of  the  month  at  the  beginning  of  the 
period  of  consideration. 

Equation  (b^)  mus t  be  solved  for  "x",  which  defines  the  period 
under  consideration.   Where  the  period  under  consideration  spans 
two  consecutive  years,  the  numerical  values  assigned  to  the  months 
in  the  succeeding  year  should  be  assigned  a  value  of  12  plus  their 
normal  value;  that  is,  January,  m  =  13-   The  quality  of  water  in 
the  reservoir,  and  consequently  the  quality  of  the  outflow,  may  be 
assumed  equal  to  the  weighted  average  of  the  measured  (or 
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calculated)  quality  of  the  inflows  from  the  tributaries  during  the 
period  of  consideration.  The  following  equation  is  developed  to 
calculate  the  specific  conductance  according  to  Method  2,  the 
running-average  method,  based  on  outflow. 


^    XZ  '=^'ji 

J  =  1   i  =  X 
C(2)„  =   Eq.  (B5) 

n       m 

0  =  1   i  =  X 

where: 

C(2)  =  EC  determined  by  Method  2; 

c   =  monthly  EC  of  a  tributary,  taken  from  the 
C/Q  equation; 

q  =  mean  monthly  flow  of  the  tributary  in  cfs; 

E  =  average  monthly  evaporation  rate  in  cfs; 

n  =  number  of  tributaries; 

m  =  calendar  number  of  the  month  under  consideration; 

X  =  calendar  number  of  the  month  at  the  beginning  of 
the  period  of  consideration  from  Equation  (b5); 

i  =  any  month  during  the  period  of  consideration;  and 
j  =  any  tributary. 
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APPENDIX  C 
HISTORIC  AND  FUTURE  LAND  USE 


The  current  water  quality  model  discussed  in  this  report  depends 
for  its  operation  on  information  developed  through  several  side 
models.   One  of  these  is  the  CU-2  model  that  provides  basic  in- 
formation on  the  consumptive  impairments  of  various  Sacramento 
Valley  stream  systems  that  are  created  by  changes  in  agricultural 
and  urban  land  use.   Land  use  data  presented  in  the  following 
tables  are  those  that  were  used  as  input  to  the  CU-2  model. 

Land  use  information  for  water  year  1955  and  I961  were  developed 
by  the  Department,  while  that  for  I93I,  1936,  and  1938  were  syn- 
thesized by  the  Bureau  of  Reclamation,  mainly  from  older  census 
data.   The  future  land  use  projections  to  the  year  1990  were  ex- 
cerpted from  file  data  prepared  by  the  Department's  Coordinated 
Statewide  Planning  Program.  ■ 

The  Department's  1955  and  I96I  land  use  surveys  were  made  by 
field  personnel  who  mapped  on  1:20,000  scale  aerial  photographs. 
The  data  were  later  transferred  to  1:24,000  scale  USGS  quadrangle 
sheets  where  planimetric  measurements  of  over  70  different  agri- 
cultural crops  and  numerous  related  urban  land  uses  were  made. 

Area  numbers  shown  in  the  following  tables  refer  to  depletion 
study  hydrographic  units  shown  on  Plate  1.  Areas  I5W  and  15E 
refer  to  those  portions  of  Area  15  lying  west  and  east,  respec- 
tively, of  the  Sacramento  River. 


^9. 


1931  Land  Use  Input 
(In  lOOOs  of  Acres) 


Crop 

Area  4 

Area  10 

Area  12+15W 

Area  13+15E 

Area  21 

Pasture 

20.9 

19.7 

21.6 

17.1 

29.2 

Alfalfa 

1.5 

9.7 

8.6 

12.1 

12.2 

Sugar  Beets 

-- 

— 

— 

— 

1.8 

Field  Crops 

0.9 

2.4 

13.8 

40.4 

13.1 

Rice 

-- 

~ 

50.8 

51.2 

8.0 

Misc.  Truck 

0.3 

0.5 

5.9 

4.6 

3.4 

Tomatoes 

-- 

-- 

— 

— 

2.3 

Orchard 

1.4 

14.7 

13.6 

33.2 

38.8 

Subtotal 

25.0 

47.0 

114.3 

158.6 

108.8 

Urban 

2.7 

3.0 

1.0 

5.0 

24.5 

Native 

107.3 

756.4 

953.8 

667.7 

716.0 

TOTAL 

135.0 

806.4 

1,069.1 

831.3 

849.3 

1936  Land  Use  Input 
(In  1000s  of  Acres) 


Crop 


Area  4 


Area   10       Area   12+15W       Area   I3-H5E       Area   21 


Pasture 

20.3 

19.3 

25.0 

19.1 

25.0 

Alfalfa 

1.5 

9.5 

9.9 

14.0 

10.5 

Sugar  Beets 

-- 

— 

— 

-- 

1.6 

Field  Crops 

0.9 

2.3 

12.9 

32.8 

11.3 

Rice 

-- 

— 

45.8 

46.2 

7.0 

Misc.  Truck 

0.3 

0.5 

5.8 

4.8 

3.0 

Tomatoes 

— 

— 

— 

-- 

2.0 

Orchard 

1.2 

14.4 

14.6 

38.1 

33.4 

Subtotal 

24.2 

46.0 

114.0 

155.0 

93.8 

Urban 

2.7 

3.0 

1.0 

6.0 

27.3 

Native 

108.1 

757.4 

954.1 

670.3 

728.2 

TOTAL 


135.0 


806.4 


1,069.1 


831.3 


849.3 
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1938  Land  Use  Input 
(In  1000s  of  Acres) 


Crop 

Area  4 

Area  10 

Area  12+15W 

Area  I3+I5E 

Area  21 

Pasture 

19.7 

19.3 

26.8 

19.4 

23.6 

Alfalfa 

1.4 

9.5 

10.5 

14.2 

9.9 

Sugar  Beets 

— 

— 

-- 

— 

1.5 

Field  Crops 

0.9 

2.3 

13.1   - 

32.0 

10.6 

Rice 

-^ 

— 

41.2 

44.8 

6.0 

Misc.  Truck 

0.3 

0.5 

5.9 

4.9 

2.8 

Tomatoes 

-- 

— 

— 

— 

1.9 

Orchard 

1.2 

14.4 

15.4 

38.8 

31.5 

Subtotal 

23.5 

46.0 

112.9 

154.1 

87.8 

Urban 

2.7 

3.0 

1.0 

6.0 

28.7 

Native 

108.8 

757.4 

955.2 

671.2 

732.8 

TOTAL 

135.0 

806.4 

1,069.1 

831.3 

849.3 

1955  Land  Use  Input 
(In  1000s  of  Acres) 


Crop 

Area  4 

Area  10 

Area  12+I5W 

Area  13+15E 

Area  21 

Pasture 

21.0 

41.2 

52.5 

43.8 

54.5 

Alfalfa 

1.5 

17.2 

43.4 

15.8 

18.7 

Sugar  Beets 

— 

0.9 

11.4 

6.5 

1.5 

Field  Crops 

0.8 

21.0 

13.5 

58.4 

14.9 

Rice 

— 

0.3 

120.8 

129.7 

49.2 

Misc.  Truck 

0.2 

0.6 

2.7 

7.9 

6.7 

Tomatoes 

-- 

— 

4.1 

5.6 

2.3 

Orchard 

1.1 

25.9 

18.2 

73.9 

43.6 

Subtotal 

24.6 

107.1 

266,6 

341.6 

191.4 

Urban 

8.4 

9.4 

2.8 

14.2 

79.2 

Native 

102.0 

689.9 

799.7 

475.5 

770.1 

TOTAL 

135.0 

806.4 

1,069.1 

831.3 

849.3 
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1961  Land  Use  Input 
(In  1000s  of  Acres) 


Crop 

Area  4 

Area  10 

Area  12+15W 

Area  13+15E 

Area  21 

Pasture 

22.3 

54.2 

54.7 

29.1 

62.9 

Alfalfa 

1.5 

16.2 

38.1 

16.2 

16.5 

Sugar  Beets 

-- 

— 

14.8 

8.0 

1.8 

Field  Crops 

0.9 

20.7 

25.0 

6.3 

13.7 

Rice 

— 

— 

111.9 

113.2 

38.3 

Misc.  Truck 

0.3 

— 

3.0 

3.9 

3.8 

Tomatoes 

-- 

-- 

6.3 

9.5 

2.6 

Orchard 

1.8 

32.8 

26.5 

84.8 

39.9 

Subtotal 

26.8 

123.9 

280.3 

271.0 

179.5 

Urban 

11.6 

11.4 

6.0 

15.8 

104.0 

Native 

96.6 

671.1 

782.8 

544.5 

565.8 

TOTAL 

135.0 

806.4 

1,069.1 

831.3 

849.3 

1990  Land  Use  Input 
(In  1000s  of  Acres) 


Crop 


Area  4 


Area   10       Area   12+15W       Area   I3+I5E 


Area  21 


Pasture 

32.0 

47.8 

37.1 

24.3 

68.2 

Alfalfa 

3.1 

31.2 

36.8 

19.8 

14.8 

Sugar  Beets 

-- 

5.3 

20.4 

11.1 

2.8 

Field  Crops 

2.4 

33.7 

52.8 

45.6 

19.7 

Rice 

-- 

-- 

179.3 

192.5 

49.3 

Misc.  Truck: 

0.7 

1.2 

6.1 

6.8 

4.7 

Tomatoes 

0 

1.4 

8.3 

14.8 

4.1 

Orchard 

5.3 

56.9 

43.0 

121.3 

64.5 

Subtotal 

43.5 

177.5 

383.8 

436.2 

228.1 

Urban 

28.1 

19.7 

6.0 

23.9 

162.2 

Native 

63.4 

609.2 

679.3 

371.2 

459.0 

TOTAL 


135.0 


806.4 


1,069.1 


831.3 


849.3 


72 


THIS  BOOK  IS  DUE  ON  THE  LAST  DATE 
STAMPED   BELOW 


BOOKS  REQUESTED   BY   ANOTHER   BORROWER 
ARE  SUBJECT  TO   RECALL  AFTER  ONE  WEEK. 
RENEWED   BOOKS  ARE  SUBJECT  TO 
IMMEDIATE  RECALL 


WB   3    'W 


LIBRARY,  UNIVERSITY  OF  CALIFORNIA,  DAVIS 

D4613   (12/76) 


3  1175  00574  5107 


L 


